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I. 


INTRODUCTION 


This report summarizes the principal findings of the theoretical and 
data -interpretation work performed under NASA Headquarters contract NASW- 
3327 during the performance period January 1980 to December 1980. The 
main body of this report is primarily devoted to an overview of the contract 
work in the various electrodynamical areas of upper and lower atmosphere 
coupling. Scientific papers published, in press or abstracts of papers in 
preparation are listed as exhibits in the latter part of this report. The 
readers interested in detailed scientific findings are referred to the exhibits. 
Some of the exhibits refer to papers in the final stages of preparation as this 
report is being prepared. For such papers, full or partial support under this 
contract will be explicitly ackxiowledged. 

Work performed under the subject contract in the following four areas 
of electrodynamical studies in upper and lower atmosphere coupling are: 

1. Magnetosphere-ionosphere-atmosphere coupling in auroral 
electrodynamics (3 papers). 

2. Middle atmosphere electrodynamics (1 paper). 

3. Thermosphere-troposphere coupling (1 paper). 

4. Tropospheric electrodynamics (progress report). 

Observational and theoretical understanding of the near-earth space 
environment have progressed in recent years to the point where it is now 
possible and desirable to find out the interrelationships between various 
components of the earth's atmosphere. In short, it is increasingly evident 
that the division of the earth's atmospheric environment into sharply defined 
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non-interacting spheres has served its earlier function of conceptual simpli- 
fication and must now be reexamined in regard to the earth's atmospheric 
response to solar stimuli. 

The major tasks detailed for the present study are primarily investi- 
gations -in-depth of studies which seem to be most fruitful in the exploratory 
work of a previous contract, NASW-3120. Studies of ionosphere-magnetosphere 
coupling in the aurora is a prime example. However, new and unique areas of 
observations bearing on the electrodynamics of the atmosphere have not been 
ignored. Items 2 and 4 above are such examples. In the area of auroral 
electrodynamics, investigations under the present and previous contracts 
have in some sense borne fruit. Work performed under NASW-3120 and 
NASW-3327 (the present contract) was the basis from which a full-scaled 
study of auroral electrodynamics at the Space Sciences Laboratory and UCLA 
was organized. This studv is presently performed under the auspices of 
NASA's Solar-Terrestrial Theory Program. 

II. AURORAL ELECTRODYNAMICS 

As a result of work performed under this contract and its predecessor 
NASW-3120, we were invited to give a review- talk, entitled "Effects of Auroral- 
Particle Anisotropies and Mirror Forces on High-Latitude Electric Fields" 
(Chiu, Cornwall and Schulz), at the 1980 Chapman Conference on Formation of 
Aurora’. Arcs in Fairbanks, Alaska. The conference is sponsored by the 
.■American Geophysical Union and the paper will be part of an AGU monograph 
(ed. S-I. .\kasofu) on the formation mechanisms of auroral arcs. 
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Auroral arcs result from the acceleration and precipitation of ring 
current and/or plasma-sheet plasma in narrov. regions characterized by 
strong electric fields both perpendicular and parallel to the earth's magnetic 
field. The various mechanisms that have been proposed for the origin of 
such strong electric fields are not mutually exclusive. However, for most 
proposed mechanisms (with the possible exception of double layers of Debye 
length parallel scale), the effects of auroral particle anisotropy, of mirror 
forces due to the inhomogeneous geomagnetic field, of auroral electron back- 
scatter by the atmosphere, and of electron trapping by the combination of 
magnetic mirroring and electrostatic forces must be taken into account in 
simulations of auroral electric fields. In addition, the effects of the very 
strong perpendicular electric field must also be taken into account in a 
kinetic description of the Poisson equation in order to achieve a unified 
theory of the auroral electrostatic structure. In this paper, progress in 
these areas in the past few years will be reviewed. It is shown that particle 
anisotropies and mirror forces can account for most electrostatic features 
of the quiet arc, while additional effects may be taking place in strong events 
in which the parallel potential drop exceeds —10 kV. 

A basic problem of all mechanisms of electrodynamic coupling 
between upper and lower atmospheres concerns how electric fields are 
transnnitted between the magnetosphere and the atmosphere. In the technical 
jargon of the field, this is known as electric field "mapping". If magnetic 
field lines were also electric equipotentials (i. e. , no potential drop along the 
magnetic field line), the "mapping" of electric fields is a simple geometrical 
procedure dictated by the magnetic field configuration only. However, as 
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can be seen from the review summary above, there exists a parallel potential 
drop in the auroral region; thus, the mapping problem takes on entirely differ- 
ent forms. We have made an in-depth study of this problem in a paper 
entitled "On the Structure of Auroral Acceleration Potentials" (Chiu, Newman 
and Cornwall). At the time of preparing this final report, this paper is going 
through the final stages of preparation; therefore, we shall include only the 
title and summary abstract in this report. Support from NASW-33 27 will be 
acknowledged. 

We examine the self-consistency of magnetospheric and ionospheric 
boundary conditions of a kinetic model of magnetospheric-ionospheric electro- 
dynamical coupling in the aurora formulated by Chiu and Cornwall (1980). 

This model includes the kinetic current conservation model of Lyons (1980); but 
in addition, it demands self-consistency between kinetic charge distributions 
and the divergence of the electric field in accordance with Poisson's equation. 

It is found that the structures of the ionospheric electric potential in response 
to imposed magnetospheric structure takes a variety of shapes in agreement 
with the "V-shaped" and "S-shaped" structures measured by electric field 
measurements on board the S3-3 satellite. In the absence of a parallel 
potential drop, the "mapping" of electric fields from the magnetosphere to 
the ionosphere is strictly dictated by the geometry of the magnetic field. 

With the presence of a kinetic model parallel potential drop, the "mapping" 
of electric fields between the magnetosphere and ionosphere is formulated 
here in terms of our Green's function treatment. It is show-n that the new 
"mapping" filters out structure scales larger than an inverted-V scale (50- 
150 km) determined by the ionospheric Pedersen conductivity. All structure 


4 


scales smaller than this natural scale can be ’’mapped through" from the 
magnetosphere to the ionosphere. 

Aside from the above two papers, we have also completed the major 
part of the scientific work on a theory of auroral arc pulsation (Chiu and 
Cornwall). We show that ion heating by electromagnetic ion cyclotron waves 
partially destroys the auroral acceleration electric, potential, giMiij; rise to a 
relaxation oscillation of the auroral potential structure. At the time of pre- 
paration of this final report, the scientific calculations have been completed 
but the paper is not yet completely written. Partial support under this con- 
tract will be acknowledged. 

III. MIDDLE ATMOSPHERE ELECTRODYNAMICS 

A major problem in middle atmosphere electrodynamics is how 
thunderstorm electric fields propagate. For radio frequency electric fields, 
the problem is well-understood. The problem for electrostatic fields has 
also beer investigated fairly completely in the last few years, due to interest 
in electrodynamic upper and lower atmosphere coupling. However, because 
the natural electrodynamic time constant in the middle atmosphere is a few 
seconds, the significant phenomenon in middle atmosphere electrodynamics 
is the propagation of sferic^ of a few second time constant. Data in this 
area have only recently been obtained with electric field measurements on 
stratospheric balloons. 

As a result of our observational and theoretical efforts, we have 
been invited to write a chapter, entitled "Sferics in the Stratosphere" (Holz- 
worth and Chiu), in the Handbook of Atmospherics, edited by H. Volland. 
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This book is under production at the time of preparation of this report. The 
theoretically predicted features of middle atmospheric sferics are in good 
agreement with observations. The essential point of this study is that inter- 
pretation of middle atmospheric electric field data must take into account 
propagation effects because the time constant for these effects can be quite 
long compared to the time constanvs of possible sources of such electric 
fields, whether they are of tropospheric or ionospheric origin. 

IV. THERMOSPHERE- TROPOSPHERE COUPLING 

We have continued compilation of a comprehensive data set of thermo- 
spheric density measurements from cold-cathode ion gauges on Atmospheric 
Explorer and USAF satellites. Fluctuations ox these density profiles in the 
thermosphere are assumed a priori to be signatures of thermospheric gravity 
waves. We have been interested in correlating features of this data set to 
physically plausible features of the troposphere as a means of studying upper 
and lower atmosphere coupling by wave propagation. Under NASW-31Z0, a 
paper was published on correlation of thermospheric wave structures with 
seasonally intense regions of deep tropospheric convective activity indicated 
by satellite-observed lightning frequency. In our present study, a paper is 
to be published (essentially accepted) on a correlation of thermospheric wave 
occurrence with transient tropospheric activity indicated by the vorticity 
area index (VAI). 

Based on the occurrence frequency of thermospheric wave structures 
in some 24, 000 measurements of density, we demonstrate that there is no 
significant hemispherical correlation between thermospheric waves and the 
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tropospheric vorticity area index in the epoch 1974-1976. Further, wave 
occurrence frequency does not exhibit the characteristic feature of the VAX 
in relation to interplanetary magnetic sector boundary crossings. In the 
polar zone, (60° - 90°' N, wave occurrence frequency does increase signi- 
ficantly if the hemispherical VAX increases above '•30x10^ km^. In con- 
trast to previous work (Rice and Sharp, 1977; Chiu et al. , 1979), which 
correlated thermospheric wave occurrence to geographically and seasonally 
persistent features of the troposphere, this work focuses upon the study of 
transient (weather) systems as represented by the VAX. Consequently, the 
null relationship between VAX and the wave occurrence frequency outside the 
polar zone may reflect the dominance of the persistent sources in the non- 
polar troposphere. 

V. TROPOSPHERIC ELECTRODYNAMICS 

Due to delays in the delivery of g).obal lightning sensor data from 
Aerospace instruments on board two DMSP satellites, we have not been 
able to complete acquisition of enough statistics until late in the year to 
correlate global lighming occurrence with solar and galactic cosmic ray 
events. Exploration for such possible connections with high statistics data 
is prompted by recent speculation that mesospheric and tropospheric electro- 
dynamics may have some relationship to nucleation around increases of seed 
ionization or to the global electrica’. circ«,ut. 

By September 1980, however, data delivery difficulties have been 
ironed out; and, we have been able to analyze the following set of lightning 
data which involve interesting solar-galactic cosmic ray events and global 
lightning events: 
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Dat?. Periods: Sep. 6-30, 1977 

Nov. 7-12, 1977 
22-30, 1977 
Feb. 13-24, 1978 
Apr. 6-26, 1978 
Nov. 2-30, 1978 
Apr. 3-30, 1979 
May 1-7, 1979 
15-28, 1979 
Aug. 14-28, 1979 

The selection criteria of these data periods aru: 1) existence oi cosmic 
ray events during the periods, and 2) existence of essentially continuous 
global lightning data from our instrxunents. 

At the time of preparing this final report, we have completed analysis 
of the above data set. A quick-look analysis shows that indeed there are some 
distinct coincident events that are mere than lo over the statistical average, 
but the total picture seems to be much more complex. We shall make a 
careful assessment of the statistical significance of the quick-look results. 
NASW-3327 will be acknowledged in our publications as a mainstay of support 
in this analysis. We plan to present the results of the analysis in the spring 
of 1981. 

VI. CONCLUSIONS 

Our research work, in a year's ti: ae, has spanned virtually the 
entire spectrum of currently interesting areas of electrodynamical studies 
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in upper and lower atmosphere coupling. Clearly, our findi'^gs cannot be 
said to have proved the existence of an incontrovertible coupling mechanism 
between the magnetosphere and troposphere. However, the work on magne- 
tosphere-ionosphere coupling, when compared with observations, indic^.ces 
that, indeed, electrodynamics are the key mechanism in which the ionosphere 
actively participates in the precipitation of solar wind energy into the high- 
latitude atmosphere. There, something mysterious happens, so that occur- 
rence statistics of thermospheric turbulence in the form of gravity waves at 
high latitudes is somehow correlated with tropospheric transient turbulence 
as indicated by ♦he vorticity area index. We cannot as yet determine whether 
this correlation is entirely due to tropospheric driving or is "helped along" 
by the special electrodynamical conditions set up in the high-latitude region. 
In this regard, it is of interest to note that the high-latitude zone (>60 lati- 
tude) is mostly outside of the most violent zone of trop.^apheric transient 
turbulence, the front of "polar" air masses. 

To find out t e possible physical mechanisms, which may be respon- 
sible for such "mysterious" correlations, we feel that the most potentially 
rewarding area of new research contributions may be in middle atmosphere 
electrodynamics. The theoretical concepts and methods in this area have 
yet to be tested by electrodynamical measurements on stratospheric balloons 
and, in the near future, on mesospheric satellites. Our contract work in 
this area has demonstrated that some element of theoretical interpretation 
of such data can be understood. We intend to pursue this area of research 
in co-ordination with balloon and satellite programs of mesospheric and 
stratospheric electrodynamics. 


The scientific papers (published, in press, accepted for publication 
or in preparation), which are attributed to support by NASW-3327, are listed 
in the following. The actual papers, presently available, are attached at the 
end of this report as exhibits. 

1. Effects of auroral«particle aniosotropies and mirror forces 
on high-latitude elec^r ■ fields, Y. T, Chiu, J. M. Cornwall, 
and M. Schulz, Invited paper, 1980 Chapman Conference on 
Formation of Auroral Arcs, Fairbanks, Alaska. To appear in 
AGU monograph 1981, ed. S. -L Akasofu. 

2. Sferics in the stratosphere, R. H. Holzworth and Y. T. Chiu, 
Chapter 16, Handbook of Atmospherics , ed. H. Volland, in 
production, 1980. 

3. Correlative study of thermospheric gravity waves and tropospheric 
vorticity area index, Y. T. Chiu and L. R. Sharp, Geophysical 
Research Letters, £, 281, 1980. 

4. On the structure of auroral acceleration potentials, Y. T. Chiu, 

A. L. Newman and J. M. Cornwall, in final stages of preparation 
and to be submitted to J. Geophys. Res. , 1981. 

5. Theory of auroral arc pulsation, J. M. Cornwall and Y. T. Chiu, 
in preparation and to be submitted to J. Geophys. Res . , 1981. 
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EFFECTS OF AURORAL-PARTICLE ANISOTROPIES AND 


MIRROR FORCES ON HIGH-LATITUDE ELECTRIC FIELDS^ 

Y. T. Chlu» J* M* Cornwall**, and Michael Schulz 

Space Sciences Laboratory 
The Aerospace Corporation 
P. 0. Box 92957 
Los Angeles, CA 90009 

Permanent Address: Department of Physics, 

University of California, Los Angeles, CA 90024 

Abstract 

Auroral arcs result from the acceleration and precipitation of ring 
current and/or plasma-sheet plasma In narrow regions characterized by strong 
electric fields both perpendicular and parallel to the earth' s magnetic 
field. The various mechanisms that have been proposed for the origin of such 
strong electric fields are not mutually exclusive. However, for most pro- 
posed mechanisms (with the possible exception of double layers of Debye length 
parallel scale), the effects of auroral particle anisotropy, of mirror forces 
due to the Inhomogeneous geomagnetic field, of auroral electron backscatter by 
the atmosphere, and of electron trapping by the combination of magnetic 
mirroring and electrostatic forces must be taken Into accoun i simulations 
of auroral electric fields. In addition, the effects of the very strong 
perpendicular electric field must also be taken into account in a kinetic 
description of the Poisson equation in order to achieve a unified theory of 
Che auroral electrostatic structure. In this paper, progress in these areas 
in Che past few years will be reviewed. It is shown that particle anistro- 
pies and mirror forces can account for most electrostatic features of the 

quiet arc, while additional effects may be taking place in strong events in 

which the parallel potential drop exceeds - 10 kV . 

*lnvited Paper; 1980 Chapman Conference on Formation of Auroral Arcs, 21-24 

July 1980, Fairbanks, Alaska 


This work is supported by NASA under Grant NASW-3327. 
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I. INTRODUCTION 


The latter half of the nineteen seventie?) has witnessed a spectacular 
increase in the understanding of the electrodynamics of the auroral arc. This 
happy circxjmstance is supported on the one hand by the advent of simultaneous 
high-resolution particle and field observations at altitudes of ^ IR^ , such 
as chose obtained by Instruments on board the S3-3 satellite, and on the other 
hand by earlier rocket and radar backscatter observations at ionospheric alti- 
tudes. Although the optical aurora has stirred human imagination since time 
immemorial, it is perhaps not an exaggeration to say that only now are human 
beings able to grasp an outline of the mechanisms for the formation of the 
auroral arcs thi‘'”"“''”e3. As we may reasonably expect, the observations 
indicate a complex but basically electrodynamical interaction between plasma 
particles and fields (electric as well as magnetic) in the auroral magneto- 
sphere. 

In this paper, we shall undertake to review certain unavoidable, indeed 
"inevitable", aspects of the electrodynamical interaction between particles 
and fields in the auroral region of the magnetosphere. These effects, which 
conspire to support an electric potential drop along auroral field lines to 
accelerate magnetospheric electrons downward and ionospheric ions upward, are 
"inevitable" in the sense that they are imposed by the following unavoidable 
external ci rcumstances: the mirroring motion of particles in the earth's 
dipole-like magnetic field, the differential anisotropy of Injected plasma- 
sheet electrons and ioas, the existence of the ionosphere with the inevitable 
implication of electron backscatter and current closure, and, last but not 
least, the existence of large electric potential drops perpendicular to the 
magnetic field. Ue shall show that an integrated model of these "inevitable" 
factors can account for the recent satellite observations of particles and 
fields in the quiet auroral arc (observed on the day-to-day basis). It is not 
our intention in this review to treat the electrodynamics of complex specta- 
cular auroral events in which electron energies are greater than a few 
kilovolts. 

Although the specific purpose of this brief review is to address the 
effects of the above-mentioned factors in auroral dynamics, it would be un- 
physicai to discuss these factors without considering their relationship to 
observations and to the theoretical interpretation of those obseirvations. In 
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particular, we need to clarify the relationship of these "inevitable" effects 
to present electrodynamical theories of the formation of auroral arcs; 
thereby, we hope to formulate a "unified" model of auroral electrodynamics « 
For these reasons, brief reviews of observations and data interpretation will 
be included here. 

In a brief review such as this one, it is inappropriate to give an 
exhaustive reference list; therefore, we shall attempt to give only a repre- 
sentative reference list, emphasizing the latest results, so that the reader 
can trace a more complete reference if he so chooses. Thus, works referenced 
here are not necessarily judged by us to be of special significance; 
similarly, works not referenced here are not judged to be comparatively any 
less significant. 


II. OBSERVATIONS 

Before the advent of high-resolution observations by auroral satellites, 
data on ionospheric and magnetospheric coupling depended on balloon, rocket, 
and radar observations which were necessarily episodic; nevertheless, the 
basic physical properties of auroral ionospheric currents and electric fields 
(e.g., Cloutier, 1971; Mozer and Manka, 1971; Vondrak et al. , 1971), together 
with their relationship to high- latitude convection electric fields (e.g., 
Cauffman and Gumett, 1972; Heppner, 1972) measured by satellites, have been 
established. Generally, these measurements have indicated that the substorm 
convection electric field in the magnetosphere drives perpendicular 
ionospheric currents consistent with ionospheric perpendicular electric fields 
(meridional and zonal) of the order of tens of millivolts/meter. In addition, 
Birkeland currents parallel to auroral magnetic field lines seem to have been 
observed (e.g., Armstrong and Zmuda, 1973), These low-altitude observations 
are primarily concerned with the morphology and large-scale processes of 
auroral substorms and are instrumental in emphasizing the importance of the 
electric field in auroral processes. They, however, have relatively little to 
say about the microscopic processes taking place in the auroral region. 

Even as the large-scale auroral processes were being unraveled, certain 
microscopic features of auroral low-energy particle precipitations were being 
discovered. Frank and Ackerson (1971) noted chat occasionally observations of 
low-energy (tens of keV) electron precipitation would show an "inverted-V" 


scruccure on an energy-time spectrum plot, i«e*, the precipitating electron 
energy spectrum hardens and then softens as the Injun 5 satellite moves 
through the structure. Evans (1974; 1975) convincingly demonstrated that 
rocket measurements of auroral low-energy electrons Indicated downward moving 
electron beams at keV energies, comparable to those of "Inverted-V* struc- 
tures. Further, by a careful study of electron backscatter from the atmos- 
phere, Evans demonstrated that these auroral electron beams are indications of 
electric potential drops, along the magnetic field, existing between the 
equator and the ionosphere. At about the same time, observations of singly 
ionized energetic O'*' ions in the magnetosphere (Shelley et al., 1972; Sharp 
et al., 1974) also gave indication that microscopic processes in the aurora 
couple the ionosphere with the magnetosphere* 

These observations of ‘‘invertad-V' structures, of electron beams, and 
of 0*^ ions of probable ionospheric origin in the magnetosphere presage very 
interesting microscopic processes to be discovered in the auroral process in 
which Che ionosphere plays an active rather than passive role. However, 
because of the episodic nature of rocket observations and because of the low 
resolution and low data rate of Che early satellite observations, the scope of 
and inter-relationship between these phenomena were not understood until Che 
launch of the polar-orbiting auroral satellite S3-3, which intercepts auroral 
field lines at altitudes up to — 3000 km , precisely in the region where iono- 
spheric and magnetospheric plasmas are expected to interact. Included in the 
S3-3 payload are instruments Co uieasure electric fields (Mozer et al., 1977), 
low energy electrons (Mizera et al., 1976), energetic ions (Shelley et al., 
1976), and plasma waves (Kintner et al., 1973). 

The S3-3 observations not only confirmed Evans' observations of downward- 
moving field-aligned electron beam at keV energies but also revealed the exisr 
tence of upward-moving ion beans aligned with the magnetic field in "inverted- 
V** structures. This certainly indicates Chat the three phenomena are 
intimately related, but more importantly the 33-3 particle observations leave 
little doubt that an electric potential drop of several to tens of kilovolts, 
aligned with the magnetic field, exists between the ionosphere and the 
magnetospheric equator. Electrostatic field measurements also indicate paired 
regions of oppositely directed perpendicular electric fields, with latitudinal 
scale lengths of some 50 km, reflecting a negative space-charge region pre- 


sumably associated with downward**screaiolag electrons* Figure 1 which Is a 
composite of particle and electric field data Illustrates the above points* 
For further emphasis, Figure 2 shows an enlarged view of the perpendicular 
electric field data for the time period marked by the brace In the middle of 
Figure 1* A cniclal, but seldom emphasized, feature which is brought out by 
the high sensitivity and high resolution of the measurements Is that the 
above correlated features are observed at the auroral zone pass after pass at 
all satellite altitudes* In other words, these are fundamental features of 
the aurora rather than episodic curiosities* 

III. THEORETICAL INTERPRETATIONS 

The key theoretical Issue concerning the Interpretation of the S3~3 
observations really Involves the electrodynamics of the auroral arc Itself* A 
key fact which must be recognized Is that a magnetlc~fleld-aligned electro- 
static potential difference of kilovolt magnitude exists between the Iono- 
sphere and the equator* In many respects, this feature has been anticipated 
In a number of theoretical considerations based on earlier observations. The 
crucial question is whether the mechanism for the buildup of such a field- 
aligned potential drop Involves the Ionosphere, for there is no doubt that the 
energy source of the aurora is derived from the Interaction between the solar 
wind and the magnetosphere via plasmas which are injected by substorm 
processes onto the auroral field lines. Some aspects of the theories of 
auroral field-aligned potential drop have been reviewed (Shawhan et al. , 1978; 
Hudson and Mozer, 1978), but our discussions will be primarily concerned with 
the particle anisotropy and particle mirroring aspects of these theories. 

Theories of auroral processes involving magnetic-field-aligned electro- 
static potential differences can be roughly classified into five categories, 
although they are not mutually exclusive. These are: 1) Debye-length double 
layer, 2) obliqtie electrostatic shock, 3) anomalous resistivity, 4) down- 
ward mapping of convection electric field discontinuities, and 5) magnetic 
mirroring effects of differential pitch-angle anisotropy between ions and 
electrons. These categories invoke theoretical arguments of varying degrees 
of sophistication and believability to show that kilovolt electrostatic 
potential drops may be produced in various assumed plasma distributions. Tne 
Debye-length double-layer model is sharply differentiated from the others by 









(Courtesy F. S* Mozer end R. 3. Torbert 



1C3 prediction of the scale length with which the total field-aligned 
potential difference is distributed, i.e. , the magnitude of the parallel 
electric field* 

Our brief discussion begins with the traditional view of these mech- 
anisms In isolation; than we shall proceed to consider their relationship to 
the effects of particle dynamics mentioned in the Introduction* 

The Debye-length double layer (Block, 1975; Shawhan et al*, 1978) is a 
boundary layer between unmagnetized cold plasma on one side and hoc plasma on 
the ocher* The potential drop across the layer is alleged to be — kl^/|e| and 
the layer thickness is of the order of several Debye lengths ( - 10 km); thus, 
the parallel electric field in double layers must be - (0*1 - 1) 7/m* A 
current-driven instability is usually invoked as the formation mechanism of 
double layers, which requires a field-aligned current greater Chan a certain 
threshold value* If potential drops inferred by electron beam observations at 
S3-3 altitudes as high as 8000 km are all due to double layers above the 
satellite, Chen evidently the ionosphere does not seem to be a factor in 
double layer formation* Freqxiently, based on observations of both electron 
ann Ion beams on S3-3, one may infer that potential drops exist both above and 
below Che satellite (Hizera and Fennell, 1977; Croley et al*, 1978). Since it 
is highly improbable chat the satellite Jtxst happened to pass through within 
the double layer thickness, such frequent occurrences seem to require more 
than one double layer to be formed on the same field line. Further, given the 
small scale length of the layer, it would be rather difficult for the mech- 
anism to explain Che coincident observation of downward electron beams and 
upward O'*' beams unless the double layer occurs at the topside ionosphere 
where 0*^ is plentiful* Theories of double layer formation are mathemati- 
cally difficult, even for very simple plasma distributions (Montgomery and 
Joyce, 1969), and a quantitative theory has yet to be develo, «;d for auroral 
plasmas in an inhomogeneous magnetic field* For these reasons, we shall not 
discuss Che possible relationship between Debye-length double layers and the 
effects of particle anisotropies and mirror forces, although the latter are 
bound to enter iuCo consideration of Debye-length double layers in magnetized 
plasmas . 

Oblique electrostatic shocks (Swift, 1975; 1976; K-an, 1975) are similar 
to double layers except chat they recognize the influence of the magnetic 
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field and consider that the shock normal is at an arbitrary angle a to the 
magnetic field direction. For the shock thickness t is measured in units 
of the ion gyroradius of a few km, a typical cross~field scale being some 15- 
20 gyroradii ( - 100 km). The field-aligned scale length is l/cos a which can 
be quite extensive if the shock normal angle a approaches x/2 . Swift has 
shown that self-consistent oblique shock solutions to Poisson's equation can 
be obtained with simple plasma distributions in an inhomogeneous magnetic 
field (Swift, 1979). The oblique shock geometry has certain advantages over 
the current-driven double layer in regard to the Interpretation of S3-3 data, 
even chough the theory was conceived prior Co S3-3. This is because the 
field-aligned scale length i/cos o can be chosen to be of the order of 1-2 
so chat only one shock (or a pair of double reverse shocks) need be invoked to 
explain the existence of potscCial drops above and below the satellite. It 
is, of course, a disadvantage that Che theory as- developed by Swift does not 
predict u , or equivalently the cross-field scale length. As we discuss 
later, this scale length can be estimated by incorporating ionospheric phy- 
sics. An oblique shock with parad-lel scale length > 1 R^ is almost certainly 
strongly coupled to the ionosphere, which at the very least supplies impor- 
tant boundary conditions for the shock* Presently, there seems to be some 
controversy concerning the distinction between electrostatic shocks and 
** double layers" of all parallel scales (Goertz, 1979; Kan, 1980). These 
arguments do not concern us; here we adopt the terminology of the initial 
author in order to maintain impartiality. 

A third mechanism by which a magnetic-field-aligned electric potential 
drop can allegedly be generated is anomalous resistivity in the field-aligned 
direction (Hudson et al. , 1973). Such anomalous resistivity may be due to a 
large number of possible modes of AC electric-field turbulence in the auroral 
plasma (e.g., Kindel and Kennel, 1971; Papadopoulos and Coffey, 1973). Hudson 
et ai. (1978) estimated that turbulent electric fields in the electrostatic 
ion cyclotron mode with amplitudes — 50 mV/m may yield sufficient anomalous 
resistivity to generate parallel electrostatic (DC) fields of — 1 aV/n. It is 
not clear how the largely perpendicular AC fields can affect parallel elec- 
tron currents »,and their resistivity). One feature common to oblique— shock, 
models and anomalous-resistivity models is that the potential smootliiy varies 
over a scale of ~ 1 ext-^ntion in order to accomnodate potential drops 
or ~ u-lu) kilovolts. It nt-st be noted that the question of how sucn an 
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extensive region of txirbtilence can be maintained at a high level (~ SO mV/m 
AC), in Che presence of non-linear stabilizing effects such as ion heating, 
must be addressed. However, to address this question, it seems to us that Che 
effects of particle motion in an inhomogeneous magnetic field and the effects 
of ionorvheric plasma cannot be ignored. Further, to study Che generation of 
Che electrostatic turbulence, the ultimate cause must be looked for in the 
characteristics of particle distribution functions. Again, we find ourselves 
returning to the question of particle distributions in an inhomogeneous 
magnetic field. 

The fourth source of auroral electric fields that accelerate ions and 
electrons in opposite directions along Che earth's magnetic field is the mag- 
netospheric convection electric field. The convection electric field is per- 
pendicular to the magnetic field at high altitudes, but its meridional Cr,9) 
component has a theoretical dlicontinuity at or near the boundary between 
closed and open magnetic field lines. Figure 3 shows the amplitudes of the 
diurnal variation of E at ionospheric altitudes (Chiu ec al. , 1979; Lyons, 
1980). Ionospheric resistivity would partially connect electrostatic equi- 
pocentials across the discontinuity, but at coo low an altitude to account 
properly for the observed component of E parallel to 3 . However, the “Itine— 
matical resistivity'* associated with magnetic-mirror forces on a hot plasma 
may Increase the altitude at which the parallel (to B ) component of E would 
appear. The details of this latter effect, which (if it occurs) would pro- 
duce the desired distribution of E • 3 with altitude, remain to be worked out 
in detail. Lyons (1980) has considered the perpendicular distribution of the 
parallel potential drop in a model incorporating the ionosphere. The effect 
would be such as to produce an upward electric field in the ?M sector (maxi- 
mal At dusk) and a downward electric field in the AM sector (aiaximal at dawn) 
of the auroral oval. This expectation la in good agreement with the diurnal 
distribution of upgolng ion beams observed by Ghielmecti et ai. (1973). How- 
ever, Co study the discribution of E parallel to 3 , one must again face the 
question of particle distributions in an inhomogeneous magnetic field. 

For the most part, these mechanisms have been considered in isolation of 
each ocher and of the ionosphere, not because pnyslclsts believe that it 
should be so, but because it is difficult to treat the couplings. In tact, a 
correct theoretical treatment of auroral phenomena will without doubt aer^e 




••v«ral of ch«t« aochanisat with ««ch other and with lonoapherlc phyaics* It 
la unfortunate that ouch of the recent literature on auroral aechaniaaa paya 
ao little attention to coupling vith the ionoaphere; aoae exceptiona (with two 
of which the authora are connected) ezlat, chough. Aa a general rule, the 
ionoaphere couplea neighboring field linea atid ullows for predictiooa of lati** 
tudinal atruccure and acale lengtha. 1. ^ *e ia no auch coupling in the 
individual aechaniaaa aentioned above (except that oblique ahocka have an 
arbitrary etructure which croaaea field linea), so none car explain arc 
structure without going beyond the given aechanisa. The unifying thread 
running through the above consideration of the aechardsaa is, of coutse, the 
necessity for establishing Che proper self-consistency between particles and 
field distributiona in an inhomogeneous mirroring magnetic field trith proper 
ionospheric and aagnetospherlc boundary conditions. 

It ia with Che consideration of this common thread in mind that we are 
led to the effects of atiroral-particle anisotropies and mirroring forces on 
Che distributiona of electric fields in the auroral region. 

IV. particle anisotropy and mirror forces 

The effects of particle anisotropies and mirror forcea in an inhomo- 
geneous magnetic field have traditionally been considered as a aechanisa for 
generating a parallel potential drop in Cbe same sense as Che mechanisms con- 
sidered above (see, for example, Leaaire and Scherer, 1974; Shawhan et al. , 
1978; Chiu et al., 1979). However, because these effects owe their existence 
to "inevitable" factors of magnetic field geometry, ionospheric and magneto- 
spheric boundary conditions, they should properly be discussed apart from 
"dynamical" mechanisms. These effects do not act in isolation or ia opposi- 
tion to the mechanisms discussed above; rather, aa we shall see, they nay be 
the common thread which unifies the various mechanisms («cept for the Debye- 
length double layer) into a single auroral mechanism. Figure 4 shows an 
example of the difference of pitch angle distrlbutlous between ions and elec- 
trons observed by instruments on board the SCATHA satellite near the equa- 
torial regions of auroral field lines. So, differential anisotropy between 
ions and electrons is observed fact, not hypothesis. 

As far as we can determine, Alfven and Falthcmmar (1963) were the first 
CO point out that low density two-component mono-energetic plasma in an 
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inhoBogancous mirroring magnaclc flald ausc support a parallal alaccrlc fiald 
unlass cba anlsocroplas of tnargy (i.a., ratio of parallal to parpandlcular 
klnaclc anargy) aca tha same for both Iona and alactrons* Thalx rasults hava 
baan furthar davalopad by Parsson (1966), yielding a potantlal drop batvaan 
tha lonoophara and tha aqua cor. 
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where B la tha magnaclc flald, v la cha kinetic energy of the component 
and spaclaa Indicated by tha subscripts, and tha subscripts 1 and 0 denoce the 
Ionosphere and tha equator raspacclvaly. Since B^^/Bq > I for auroral field 
lines, la I A# would be much graacar chan cha auroral particle kinetic energies 
( ~ 10-*100 kaV) unless cha anisotropy difference between electrons and Ions Is 
of order B^/B^ • Ponyawln ec al. (1977) pointed out that the potential drop 
for a tvn-componenc bl'-MazwalUan plasma would be otuch less. 
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where s ■ B^/Bq I, yielding |e| A^ to be about half of the average particle 
kinetic energy (a few keV). Similar results have also been discussed by 
Whipple (1977) and Lennartsson (1977), 

While the geometrical and boundary~lmposad properties of particle dlstrl** 
buclons are the basic causes of cha parallal potential drops discussed above, 
l.e«, Ions and electrons of different pitch angle distributions mirror on the 
average at different locations of the field line thus creating a charge~ 
separation electric field, the above results are arrived at by requiring a 
parallel potential distribution to maintain a balance of positive and negative 
charges on a field line (qtiasl-^eucrallcy) , For this reason, the effects of 
particle anisotropy ana mirror forces are quite frequently intermingled with 
quasi^neutrallty ~ to the extent Chat these effects have become inseparable 
from quas 1-neu trail ty. la reality, we must be careful to distinguish quasi- 
» which is nothing more or less chan an assumed method of deriving 
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the parallel poceaclel ^ , from the more basic and "inevitable" factors of 
particle anisotropy and oilrror forces* Indeed, we shall discuss aore general 
applications of these factors In the next section. 

The application of quasl'-neutrallty to obtain solutions of Iclnetlc equi- 
librium between particles and fields has had a long history. Lemalre and 
Scherer (1971) have pioneered a series of applications of these kinetic prin- 
ciples to a number of aagnetospherlc problems Involving primarily Isotropic 
Maxwellian distributions. Applications to auroral plasmas have become the 
dominant topic In recent years (Lemalre and Scherer, 1974; Knight, 1974; 
Lennartsson, 1977). Chiu and Schulz (1978) reexamined the problem In rhe 

light of the S3-3 observations and suggested the necessity of caking Into 

% 

account other "Inevitable" factors such as electron backscatter from the 
ionosphere (eloquently noted by Evans (1974; 1975) a few years earlier) and 
electron trapping by a combination of electric and mirroring forces. A 
typical quasi-neutral solution for the equilibrium potential Is shown In 
Figure S. In general, the parallel electric field due to the factors 
discussed here Is distributed over the entire auroral field line so that the 
'parallel electric field strength Is < 1 mV/m. Further, as- has been noted In 
simple model calculations, the more complex models confirm that a potential 
drop of < 10 kV can be supported by these basic factors of particle kinetics; 
thus, particle anisotropy and mirror forces cannot account for the more spec- 
tacular auroral events In which the potential drop can be as high as 30 kV. 

A final aspect of particle anisotropy and mirror forces in auroral 
dynamics is somewhat subtle and Is not generally appreciated. By suggesting 
that the basic parallel electric field distribution Is a direct kinetic conse- 
quence of externally Imposed geometric and boundary factors, we are also 
suggesting that there Is an "Inevitable" kinetic resistivity to the auroral 
low-density plasma so that dynamical processes such as plasma waves do not 
play a major role in determining the current-field relationship ("Ohm's law") 
in Che quiet arc. This viewpoint has been expressed as the "kinetic picture" 
by Lemalre and as the "current-carrying characteristics of the low density 
auroral flux tube" by Evans. A simple but "inevitable" prediction of this 
"kinetic picture" is the relationship between field-aligned current Jj and the 
electric potential drop 44 
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whare A is darlvsble from c&« psraaacars of tha parclcla dlscrlbuclon func** 
tloas (Fridman and Lamalra, 1979; Chiu and Cornwall, 1980). Lyons ec al., 
(1979) bava esaantially demonscracad tha valid! ey of (3) wlch obsarvacions. 

V. TOWARDS A UNIFIED MODEL OF THE QUIET ARC 

In Saccion III, wa hava nocad ChaC cba tcinamacics of particla distribu- 
Cions in cba inbocoganaoua mirroring gaomagnacic field is a cotmon thread 
which links Cba various suggescad aachanlsas for cba auroral parallel electric 
field. In Saccion IV, wa hava ravlawad Che davalopoenc of this common thread “ 
specifically, cba effaces of parclcla anlsocropy and mirror forces in quasi'* 
neutral models. Since Chase effaces owe chair existence to "inevlcable** 
factors of magnacic field geomacry and boundary conditions in the ionosphere 
and magne cosphere, the next logical step In Che development of a model of 
quiet arc formation would be Co consider Cba various mechanisms (except Che 
Debye-lengch double layer) in Che light of this unifying thread. 

Ic is evident chat these physical mechanisms do not exist entirely inde- 
pendently of one another. For example, if the restrictive assumption of 
strict charge neutrality is removed in the magnetic-mirror model, one has an 
oblique electrostatic shock in a mirroring field. To the extent that no 
dissipative mechanisms such as wave— particle turbulence are included in such a 
'shock'*, the resulting electric-field structure is better described as a solu- 
tion of Poisson's equation. From another point of view the oblique shock can 
be described as some sort of zero- frequency electroqcatic ion-cyclotron (ESIC) 
mode. There surely is a great deal of ESIC turbulence connected with auroras, 
and Che physical distinction between Che oblique shocks of Swift and this tur- 
bulence is at best imprecise. Yet the merging of wave turbulence and shocks 
can lead to substantial parallel electric potential drops in the complete 
absence of anomalous resistivity. (The reader need not be reminded that tur- 
bulence is not synonymous with anomalous resistivity; in fact, it is quite 
difficult to make anomalous resistivity out of even the most turbulent waves.) 

The effects of AC electric fields and electrostatic shock structures, 
together with the very important factor of ionospheric current transformation 
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c«n all be folded Into a model which Includes the effects of particle 
anisotropy and mirroring forces* The formal development Is based on charge 
conservation (Polbson's equation, of which quasi-neutrality Is a truncated 
form) and current conservation In the Ionosphere. Chiu and Cornwall (1980) 
have considered such a model In dipolar magnetic geometry, coupled with 
ionospheric physics. In such a model Che parallel potential drop Is Inti- 
mately coupled to Che perpendicular electrostatic field structure as Indicated 
In Figure 1. 

Poisson's equation with perpendicular electric field In a dipolar 
magnetic field 




I n e 

J ^ ^ 


(4) 


becomes 


K ® if (b“ E,) - 4» J SjCj 


(5) 


where K Is Che plasma dielectric constant: 


n.M.c‘ 

1 + 4v J J 


( 6 ) 


and Che density n^ Is that which would be appropriate In the absence of the 

perpendicular electric field. Equation (5) Is an approximate generalization 

of Swift's formulation (1975) to an Inhomogeneous magnetic field and Is in 

agreement with Swift's own generalization (1979). The quasl-neutrallty 

equation of Chiu and Schulz (1978) Is simply Chat Che right-hand side of 

(5) vanish. In fact, talcing Ej^ from data and V«E^« ^^® 

perpendicular scale, one finds Chat the left-hand side of (5) yields a frac- 

-2 -3 

clonal charge separation d n/n ■ 10 -10 . Or, loolced at Che ocher way 

around, this sort of charge separation can drive perpendicular electric fields 

of 0(M*’/ t^e) • lOO's of mV/m, as well as parallel fields of 0[( ]. 

It should be noted that Che term B9(B ^E,)/3s In (3) is quite negligible; it 
2 2—8 ' 

is of order ^ /1| "10 where is the Debye length. This is why ordinary 

double layers need such small parallel scale lengths in order for there to be 
appreciable charge separation. In this generalized formulation, the exis- 
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eenc* of E| ( indep«adcnc of vhechar ch« second tera on Che lefc side of (S) is 
included or noc) is supported self-consiscencly not only by particle aniso' 
tropy and alrror forces implicit in n. but also by the ezlscence of parpen- 
dicular potential drops in the magnetosphere implicit in K 
since K > 1. Thus, a discontinuity in 7*E^ in the aiagnetosphere (Lyons, 
1980; Chiu at al. , 1979) an also be contributory to E|. 

From (S), introduces the perpendicular scale length of the 

auroral arc into the model. However, to self-consistencly decermlne , (S) 
is not enough because current and charge must also be conserved at the ionos- 
phere. In ocher worl, the scale must be such that ionospheric and magneto- 
spheric conservation iivs are both satisfied. Thus, the scale length of the 
perpendicular electrostatic field structure is related noc only to Che field- 
aligned current to the ionosphere but also to the ionospheric Pedersen conduc- 
tivity. Thus, ionosphere-magnetosphere coupling is a crucial ingredient 
determining the geometric structure as well as the energetics of the quiet 
auroral arc in such a model. 

To incorporate ionospheric current transformation In the region (120- 
2000) km altitude is a very difficult cask. At present, only very rudi- 
mentary methods are used to prescribe the current flow in the ionosphere and 
Che horizontal electric-field gradients there; they are coupled by the 
requirement chat the ionosphere must be able to balance the current flow (as 
well as to balance sources and sinks of charge). Two simple equations express 
these facts: (e.g., Atkinson, 1970; Coronlcl and Kennel, 1972) the equation 

of current continuity 


J 


I 


_3 

H. 




and of charge balance at the ionosphere 


dc 



(7) 


( 8 ) 


In these equations, x is a horizontal north'south coordinate, £ is the 

? 

height-integrated Pedersen conductivity, M is the height-integrated tono- 
spheric electron density, a is a recombination coefficient, represents 

aon-aurorai sources, and S represents sources of .harge from auroral electron 
precipitation, or from upward-going ions. 
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At pr«««nt, no Mtlsfactory solution of such s aodtl has yst baan 
obtained In the rotunt~curranc raglon, although an approximate solution In the 
central electron beam region has baan obtained (Chiu and Cornwall, 1980). A 
schematic Illustration of this model Is given In Figure 6. Figure 7 shows the 
latitudinal structure of the approximate solution which can be compared with 
the observed structures shown on Figure 2. 

V. CONCLUSIONS AND PROSPECTS 

Ue have undertaken a very brief review of the various mechanisms which 
have beun invoked to account for the observed consequence of parallel 
potential drops In the auroral arc. Our purpose Is, of course, not to address 
the details of these mechanisms separately but to emphasize that all the 
mechanisms (with perhaps the exception of Debye-length double layers) have a 
commonality underlying then: the unavoidable effects of particle anisotropy 
and mirror forces, and the necessity to Include Ionospheric physics as well as 
proper magnetospherlc boundary conditions. While specific dynamical 
mechanisms may be of varying importance In the vide range of auroral condi- 
tions, It seems to us that particle enisotroples and mirror forces must be 
taken into account for most electrostatic features of the quiet arc when the 
parallel potential drop does not exceed a few kilovolts. For more 
energetic events, additional dynamical factors such as electrostatic shock and 
perhaps anomalous resistivity may come into play. 

To formulate a unified model of arc formation,' i have pointed out that 
the basic geometrical factors of particle anisotropy and mirror forces must be 
merged into the considerations of the various ‘’mechanisms’*. In the light of 
these basic factors, consideration of the various ‘‘mechanisms'* in isolation of 
each ocher has no a priori physical Justification. Indeed, it is an unfor- 
tunate historical accident. Auroral physicists are now beginning to realize 
the necessity to unify the diverse mechanisms into a comprehensive model of 
auroral arc formation. The task is difficult because it must encompass both 
magnetospherlc and ionospheric physics, but attempts in this direction are 
already being jade (some features of these preliminary attempts are reviewed 
here) . 
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7. OiscrlbucloQs of •Icctroacaclc pocaoclal U and Che norchrsouch 
perpendicular eleccroacaclc field E, as functions of L ac 2000 ka 
altitude. ^ 
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Tha alaccrodfiualcal scruccura of auroral area la obvloualy much aora 
cooplax chan cha vary aiapla faccora dlaclllad Inco cha aodala ravlawad 
bara. Tha iaaadlaca proapacc for choaa aodala would ba co conaldar cha baalc 
alacerodynaalcal acruccura in ralaclon co apaeiflc phanoaana ooa ac a elaa. 
For axaapla, cha parclcla dlacrlbuclon funcclona in aalf~conaiaeanc aqulll* 
brlua vlch cha alaccroacacle pccanclal # naad noc ba aeabla wlch raapacc co 
varloua plaaaa wava aodaa; Charafora, avaluaclon of wava ganaraclon undar cha 
baalc unavoidabla condlclona of parclcla anlaocropy and alrror forcaa In a 
strong DC alaccrlc flald will ba an laporcanc araa of auroral chaory In cha 
naar fucura. Furthar, alchough It la oucalda of cha scopa of cha praaanc 
ravlaw, cha affacca of cha plaaaa wava AC alacCrlc flalda ara bound co faad 
back upon cha DC alaccrlc flald acruccura and cha parclcla dlscrlbuclon. Tha 
naxc logical scap of arc~aodal davalopnanc will probably ba In chla dlracclon, 
chua unifying soaa aapacca of alaecroacaclc curbulanca chaory Incc tha baalc 
faccora conaldarad hara. 
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Abstract 


We examine the self-consistency of magneto spheric and ionospheric 
boundary conditions of a kinetic model of magnetospheric -ionospheric elec- 
trodynamical coupling in the aurora formulated by Chiu and Cornwall (1980). 
This model includes the kinetic current conservation model of Lyons (1980); 
but in addition, it denuinds self-consistency between kinetic charge distri- 
butions and the divergence of the electric field in accordance with Poisson's 
equation. It is found that the structures of the ionospheric electric potential 
in response to imposed magnetospheric structure takes a variety of shapes 
in agreement with the"V-shaped" and "S-shaped" structures measured by 
electric field measurements on board the S3-3 satellite. In the absence of 
a parallel potential drop, the "mapping" of electric fields from the magneto- 
sphere to the ionosphere is strictly dictated by the geometry of the magnetic 
field. With the presence of a kinetic model parallel potential drop, the 
"mapping" of electric fields between the magnetosphere and ionosphere is 
formulated here in terms of our Green's function treatment. It is shown 
that the new "mapping" filters out structure scales larger than an inverted-V 
scale (50-150 km) determined by the ionospheric Pedersen conductivity. All 
structure scales smaller than this natural scale can be "mapped through" 
from the magnetosphere to the ionosphere. 
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I. imooucnoN 


OiililM sfarle* la Cha lowar croposphara, a alarlc Man from tha atraco^ 
■phara la ooc uaually eharaecarlzad by an loolzaclon ehannal producad by a 
▼lalbla alaetrleal dlachatfga. Tha aablanc eonducclvlty abova 20 Voi alcleuda 
la ganacally aufflclane Co carry a afarle ralatad currant pulaa without alac~ 
trlcal braakdown. Thla la a raaulc both of tha hlghar aablant alaetrleal 
conductlylty aa wall aa tha larger area available for a current path aa tha 
algnal propagataa farther away froa tha aouree In tha cloud* To data vary fav 
in altu aaaauraaanta of afarlc'^elatad alee trie curranta and flalda in tha 
atratoaphara have bean reported* 

Sfarle alaetrleal aaaauraaanta made onboard aircraft flying over thxindar* 

% 

atone in tha tropoaphara have a long hlatory (ef* Halt, 1953; Stargia at al*, 
1957; and Vonoagut ec al*, 1966)* Mora racantly, atratoapberlc af erics have 
bean dataetad with hlgir*altltuda~balloon-borna electrical instruaentation by 
Benbrook. at al* (1974), Burke (1975) and Holrworth (1980). 

ThaM three racene works have concantrated on thunderatona related strat'~ 
oapnarlc electric fields in general and each reported briefly on stratospheric 
sfarle signatures. This chapter will extend these and tha tarller works by 
describing stratospheric sfarle related electric field detection oathoda and 
present absarvatlons of a variety of different signatures* Propagation of the 
sfarle source disturbance in tha stratosphere from which a physical descrip* 
cion of Che observed claa scale and slgnacures can be detaralned, will be 
aachematieally described. 
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11. STRATOSPHERIC ELECTRIC FIELD INSTROMENTATIOM 

Tha dacaetor oaad la tha azpacloanca daseribad balov is a doubla laagBoir 
proba spaeifleallj daslgnad to oparata in cha atratoaphara (aaa Flgura 1). 
Tha davica aaaauraa Cha potaatlal diffaranca batwaan two Idanclcal. aaparatad 
eonductora «lth a dlffacantial oparatlonal aaplifiar. Tha alaetrlc fiald 
eoaponant along tha batwaan cha two probaa ia Cha pocanclal diffaranca 

dlTldad bf cba proba aaparaclon dlacanca. 

Manj chaoracleal and axparlnantal dlaeuaalona of cha doubla Langmuir 
proba (aaa Hocar, 1973, and rafarancaa Charaln) haya ahoun lea applleablUcy 
Co balloon-'boma maaauramanca ena abaoluca dc alaccrle fiald In cha acrato- 
aphara to an accuracy of a fav mllllTolea/maCar. Thla aaeclon will dlacuaa 
cha particular problcaa ancountarad uhan applying cha doublan>coba eachnlqua 
Co Cha uppar acaoaphara and daacrlba Cha phyaleal apparacua uaad In cha axpar* 
Inanca. 

Tha praaanca of amall aolacular Iona proyldaa cha acaoaphara vlch a low 
alaccrlcal conducclvley ranging on cha awaraga from 5 z 10~^^nho/a ac cha 
ground CO 10~^^aho/a ac 30 ka alclcuda. In ordar co aaka a doubla~p?oba 
pocanclal diffaranca aaaaureaanc In cha acaoaphara cha dacaccor auac hava an 
affaeclva Incamal alaccrlcal raalacanea which ia tzcraaaly high for cha 
following raaaon. Tha ^cancial diffaranca V in cha acaoaphara ttHI cauaa a 
curranc 7/R^ Co flow enrough cha alaccronlc Inpuc lapadanca 8^^^. This cur~ 
ranc recuma through Cha acaoaphara along a pach wich raalacanea R^^a* ^aca 
is coapucad froa cha raquirad pocanclal diffaranca co collacc a glyao curranc 
whara chac curranc is iapadad by cha acaoapharle raalscirlcy. Mozar and 
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FlRure 1. Stratospheric electric field detector showing three pairs of 
nutually orthogonal spherical probes. Numbers on probes are 
referenced In later figures. 



Sarlla (1969) cAleulAC* ehae for « ophorlcAl probo of radluo In eb« aeno* 

aphara with eonduetlTlCT 9, that ^ • lO^^O ae 30 taa alcicuda* 

0 

Qdlaao 1^ » eha roleaga dlTldar affaec of thaaa evo raaiatora vlll 

eauaa aoac of tba roltaga to appaar aeroaa *lta* Thua Co naaaura an unaetanu-* 
aead ^ralua of ▼ eba alaceronlca *000 aaelafy ^9t eha daeaeeor eo 

aaae eha raqulraaane of high Incarnal Inpadanea (Z > 10^^ ohas), eha inpue 
algnal fxoa aaeh saparaead eondueeor la fad Ineo a vary hi^ lopadanea FCT 
oparaclonal anpllfiar (op aap)\ In addleion eo ehia op aap ae eha Inpue, eha . 
eoaplaea Inaeruaant la daslgnad In a aannar ao aa eo radoea laakaga eur~ 
ranea* Badoeing eha laakaga la aeeoapllabad hf clreuie daalgn eooaidaeaclotio 
aa vail aa azeamal alaeerleal elaanllnaaa» Tha ahiald of eha coaxial eabla 
eo eha proba la drlvan ae eha aaaa poeaneial aa eha proba leaalf, ehacaby 
cadueing laakaga* Driving eha ahiald ae eha proba poeaneial alao allainaeaa 
capacleiva eorranca which alghe darlva from physical aociona of eha Inacru* 
aane* Tba op aapa ehaaaalvas ara aoimead oa a hi^ lapadanca eafloa basa and 
ehan alacerlcally shlaldad* Tha axeamal surfacaa of eha Inseruaanc ehac ara 
In coneace vleh eha probaa noac also bava high alaeerleal rasiseanea* Tha 
fibarglass booaa on which eha probaa ara aouncad ara carafuUy elaanad vleh 
Froon aad drlad In a claan aarlrooaane eo raaova all oils aod aolscura boch 
Inslda and out* k high Isipadanea coaxial eabla la ehan elaanad and scrung 
dovD eha eanear of eha boon and boch onda ara plugged. 

Two of eha dlffleulelas sneounearad whan asking a poeaneial dlffaranca 
naasuraaanc In eha aoosphara arise fron work funeelon dlffarancas bacasan eha 
evo probaa and fron phoeoalaeerle anlssloa enrraacs eo -cha probes. Tha work 
funeelon dlffarancas ara ainlalxad by ewo proeasaas. Tha probes ara aaelraly 
eoaead vleh a soluelon of aquaous carbon called *aquadag* which halpa co aaka 
cha work funeelon unlfom. During daea eoUacclon eha anclra apparacus la 
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roCACtd about a ▼•rtical axis. Tha diffaranct In potantlal batwean cha tvo 
probas dua to an aztarnal alactric flald vlll than appaar a# a sinusoidal 
yoltaga variation In tbs dlffaranca oaasuraoant, but any raaalnlng nat work 
function dlffaranca batvaan tha probas will than ba oaasurabla as a dc offsat 
of this sinusoidal voltaga* Tharafora, analyzing just tbs oscillatory part of 
tha output signal allows afflciant ranoval of thasa work function arrors In 
tha horizontal alactric flald naasuraaant and datarmlnas the real zero voltage 
basa line for translant cvants. 

Tha photoalactrlc affact can causa a significant current to the probe If 
tha fluz of ultravlolat (UV) radiation Is larga* At balloon float altitudes 
naar 20 kllonatars tha photo currant Is not significant because most of tha OV 
absorbing ozona la abova this altltuda. Howavar, for ezperlaants above 30 to 
35 kllonaters tha photoalactrlc currant can ba significant for saveral hours 
on althar side of the diurnal oaxlnua zenith angle. To ainitaize this problem, 
probes for high altltuda experiments are designee to have equal solar cross- 
section and are usually spherical so that the photoelectric current Is nearly 
the same for both probes. Thus the nat currant to the probes from this effect 
will produce only a «mall dc offset. 

Tha probas along the spin axis tmist make an absolute potential difference 
measurement without tha benefit of the regular payload rotation. Thus the 
vertical potential difference maasurement will be biased by work function 
offsets. This dc offset In the vertical electric field measurement can amount 
to about a ±20 mV/m uncertainty in the vertical electric field, which is small 
compared to the atmospheric electric field at 30 km altitxide and small com- 
pared to most sferlc fields. Because of this inherent vertical field uncer- 
tainty and because this bulky balloon payload is difficult to launch, it is 
sometimes desirable to remove the lower vertical boom. Then the vertical 
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«l«ccric fl«ld is dscsralnsd troa th« pocsacial dlfftrsocs bscwssn ths upp«r 
prob« and chs svsrsgt of Chs four horizoacsl probss. This oskss ehs tlgnnl Co 
aols« rsclo lo«*«r by s factor of 2 and cha probas sllghcly asynmacrle, buc cha 
lack of a lowar booa oakas Cha awkward payload in Figura 1 coosldarably aaalar 
Co launch. Tha data prasancad in this chapcar ara from payloads wlchouc Cha 
lowar varclcal proba. Tha inscruaancs flown vara coaposad of flva conduceors, 
alchar ona-fooc square flac placas or spharas aach a fooc in dlaaacar on chraa 
oucually orthogonal pairs of booms surrounding a conducting payload ground. 
Whathar spharas or cha 'lass azpansira flat places are flown depends on Cha 
aagnicuda of Cha azpaccad phocoalaccrlc curranc at cha balloon float aid'* 
Cuda. Above Cha uppacaosC varcical sphara is a small motor which is connected 
CO cha load Una and provides Cha fifteen second rocaclon period. Tha entire 
payload, including calamacry alaccronlcs, alclaacar, aspect magnacomacer, and 
power supply, weighs about 30 pounds. 

Tha potential dlffarancas bAtwaan th? three orthogonal secs of double 
probas ara calaaaterad to a ground recaivlog station by various mathods. Tha 
current payload configuration aaploys an onboard analog Co digital convartar 
from which cha subsaquant biphase encoded PCI modulation signal makes down~ 
scraam data processing simple and accuraca. 


III. OBSERVATIONS OF SFERICS OVER THUNDERSTORMS 

Much h«s b««n written About the quasl'-dc fair weather electric field in 
Che etretoaphere (cf. Moxer end Serlin. 1969; Mozer and Lucht, 1974; Bolzwortb 
et el., 1977). However, very few obaervaciona of tbunderatora related etreto"- 
apheric electric fielda have been reported. Thia aectioo will auamarize the 
electric aignaturea of aferics aa aeen in hundreda of hours of balloon data 
collected near 30 ka altitude by R. Bolzvorth and F. Mozer at the University 
of California, Berkeley. 

Figure 2 ehowe the ciaa history of a large thunderstorm in the dc elec~ 

trie field at 27 ka. This data was taken in August 1974 from a balloon flight 

out of Ft. Simpson, N.V.T., Canada. The coordinate system used for all fig~ 

urea in thia chapter has z in the eastward magnetic direction, y in the north- 

ward magnetic direction and z in the vertically upward direction. Thus a fair 

* 

weather vertical field would point toward the ground and would be negative. 
The upper panel in Figure 2 is the vertical field which shows ambient fair 
weather electrical conditions of a few hundred mv/m negative both prior to and 
following Che thunderstorm. During Che storm there appear Co be evo cells or 
two fronts passing under the oalloon with Che second larger chan the first. 
The total horizontal electric field is shown in the central panel for com- 
pleteness and cbe boCCoo panel shows the sferic race. For more detailed 
discussion of Che dc effects on all three electric field components and the 
conductivity during such a thunderstorm, the reader is referred Co Holzvorth 
(I960). In Che time history of the sferic rate one can again see the double 
peaked structure although the peaks do not coincide in time. This figure is 
to be used as a description of ch.e ambient conditions during which detailed 
individual sferlcs will be analyzed below. 
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^UG 13. 1974 


aUG 14. 1974 


^electric field over a large thunderstorm near Fort Simpson, 

• • anada (6„ N, 121.5 U). The upper panel Is the vertical 

^ tctrU leld, central panel Is the absolute value of the total 


horizontal electric Held 
sferlc rate. 


and the bottom panel Is the 


^1^ 


Fljur* 3 is Also for bsckgrouad purposss sod shows two DHSP wAschsr 

pbocogrsphs ac class which brsckac cho dscs in Flgur* 2. No oehar obsArrs- 

clonsl wAAChar dACa is ATAllAblo oehar chan dacsllad launch slCa raports of 
cba local cloud foraacloas and praclplcaclon* Tha balloon aovad ac abouc 13 
knocs so cba waachar ac cha local launch slca was noc diraccly appllcabla eo 
chac vmdar cha balloon by a faw hours afcar Launch. In any aTanc Cha aonC 

chac can ba said is chac an axcauslTa cloud syscaa did swaap frost ao.'Chwasc Co 

souchaasc acroas Cha raglon as shown In Flgura 3 and chac sfarlcs wtra ra~ 
porcad on our radio Calaaacry fraquaoclas as wall as* ground raporcs of rlslbla 
lightning and Chtndar ChroughouC auch of Cha flight. 

Flgura 1 shown a 20-nlnuta sagMnc of alaccrlc flald data during Cba 
chundarscom of Figuraa 2 and 3» Tha larga spurious signals ara dlraccly 
eorralacad eo scaclc oolsa burses ovar cha audio fraquancy celaaacry channal 

and CO eras pond co sfarlcs. Tha cop panal of Flgura 4 is cha ▼arclcal alaccrlc 

« 

flald ac loMsc gain whlla cha aacood panal froa cha cop la also cha varclcal 
flald buc now ac a gain up by 3. Tha boccoa Cwo n«aals aach show tha poccn- 
cial dlffaranca baewaan ona of two oppoolca horlzonCal probas Co payload 
ground. 74 ij ISO dagraaa away froa 73. Dua Co various gal:i channels avail*- 
able Cha horizontal ancannas ware often saturacad while cha lowesc gain verti- 
cal 'as never sacuratad in this flgura. 

Tha data fron cha hcrizoncal probas (7^ and 7^ in cha boccoa cwo panels) 
rapresanc cha individual probe voltages ralaeiva Co payload ground. They have 
bean prasanced individually in order co show cha dranaclc efface on payload 
ground as cha sfaric passes. Careful iospacCion of chase figures -w ill show 
chac oucsida of cha sfaric cruisiancs, 7^ and 7^ always have cha opposite sign 
ac any given tlaa. In tha case of every sferlc in Flgura *♦, cha signals fron*. 
cha cwo opposite horizontal probas depart from cha oasaliaa in opposite dlrtc- 
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*4. Twr.av rinuti's of r.iw t ield ilar.i i1urln\: thi» r humit'r s t orn 

Mr.tiro'*. .uul i. Tilt* two nppor panels are the vt‘rtical ele<MrU' 
t'iehl, at (wo Jitiorent >»ains. The bottom two puiels are I wt^ 

’^ivisn ronenl s ot t h.e lv>ri7ontal electric tieki from [^rt>bes ^ and 
win. h ,\re ISi'*' out o{ phase. (See Kipaire I for or i ent at i oju ^ The 
backv,round <>scillatorv electric field of the hortt^m two fip.nrt's is 
Ine to t!u* ambient Pr electric field while the spikes are <ine to 
s i cr i cs. 



Cloas Initially. Then, however, the payload ground goes alternately very high 
then very low as seen when and Vj have the same sense near the end of each 
large sferlc* This payload ground fluctuation has been removed from the upper 
two vertical panels showing the horizontal field by subtracting the average of 
all four horizontal probe signals from the upper vertical probe potential. 
This payload ground transient is seen consistently throughout our data set for 
all large sf erics but does not affect our high impedance circuitry , which is 
highly current protected. There are at least two smaller sf erics in Figure 4 
which were unaccompanied by gro^uld transients and several others in Later 
figures. The next to the last sferic in Figure 4 at about 0016 UT which did 
not cause a ground fluctuation will be evaluated more careftiUy later. 

There are many features to note from this figure. First of all, most 
sferics last for about 5 to 10 seconds and have the sense which would appear 
from the rapid raooval of positive charge below followed by the slower return 
to prestroke values. Another interesting thing to note in Figure 4 is that 
the sferics are not always related to major changes in the dc level. This can 
only be accurately determined for dc steps larger than about 100 mv/m from 
these data because of the residual rotation period fluctuation in the vertical 
data due to slight asymmetries of the verical measurement probes. However, 
small discrete dc steps are sometimes seen as shown below. 

The next three figures show other 20~minute sequences of data containing 
sferics. Figure 5 is from the same thunderstorm shown in Figures 2~4 and 
shows an epoch between the twa major activity periods of the storm. The 
interesting feature is the sferic near 0056 UT which has a positive precursor 
to the negative sferic transient. This brings up the point of temporal reso- 
lution. Ail of the data in this chapter were telemetered at a rate allowing 
complete vector field determination every 0.5 second. Thus, while the elec- 
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troaica ar® ouch faster than this, time resolution Is limited by the telamec- 
ry* It is therefore possible that more of the sferics displayed in this 
chapter actually had positive precursors but were missed because of time 
resolution* This type of precursor will be seen again in a later figure. 

Figures 6 and 7 show two succeeding lO-nninuta segments of data during a 
thunderstorm seen from a flight over Thompson, Manitoba, Canada in August 
1978. In this flight, as before, there is only a single upper vecclcal probo 
from which the average of the horizontal probes Is subtracted. However, now 
the horizontal field Is very large making the asymmetries between the payload 
ground planes and the upper vertical probe more apparent. The difference 
between and 73 In the lower panels (or Just the envelope of their signals) 
divided by probe separation gives the horizontal fields. In this case a 
horizontal potential difference of up to 8.0 volts (in Figure 6) divided by 
probe separation gives a horizontal field of over 2.5 volts/meter. This is 
over twice as large as the vertical field seen In the upper panel. Also, in 
this flgxire it is now possible to see a net dc shift of about 300 to 400 mv/m 
In the vertical field after a sferic at 2302.5 QT. The next 20-minute period 
In Figure 7 shows a lessening In severity of the storm and a gradual return to 
lower vertical and horizontal dc fields. The last two sferics in Figure 7 are 
quite Interesting In that two have an inverted polarity. The last one again 
shows a precursor but now In the negative sense compared to that in Figure 5. 

A summary of the sferic signatures seen the vertical field on these 
balloon flights Is shown In Figure 8. Of the four basic slgnatxires the most 
common In all these data Is the top one of Figure 8, namely, that which would 
be expected by the removal of positive charge below. The other three cases 
are seen only rarely with Figure 7 showing Che only example we have seen of au 
inverted polarity stroke with an opposing precursor as In the bottom signature 
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Figure 6. Twenty minutes of raw electric field data during a balloon flight 
from Thompson, Manitoba, Canada (55® 45*N, 97® 51 *W), DC shifts 
resulting from the sferics are now visible in the vertical field. 
In this case the ambient horizontal electric field is larger than 
the vertical field. 
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Igure 7, Twenty minutes of raw electric field directly following data In 

Figure 6. Note: the reversed sferlc orientation and precursor near 
the right hand end of the vertical field panel. 
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17. THEORI OF SFERIC PROPAGAIION IN THE STHAIO SPHERE 

In order co Interpret the sferlc signal structures detected by electric 
field aeesureoents In the stratosphere, it Is not only necessary to have flra 
Inforaatlon on the spatial and teaporal structures of lightning discharge 
currents, but it is also necessary to establish basic Ideas about sferlc 
propagation. The necessity to unravel these tvo factors In data Interpreta** 
tion Is simply that tlas scales of the sources are folded Into time scales of 
the propagation «"■ to produce the detected signal. While the character* 
Istica of lightning discharge currants are dealt with In Qiapter 2, the theory 
of sferic propegation in the scracoephere will be reviewed and forauJated In 
such a way thee the relatlooship between source, propagation oedlum and detec- 
ted signal will becone transparent. This entails extensive usage of the 
theory of Green's functions. 

% 

The theory of sferlc propagation in the stratosphere, over horizontal 
distances which are short coopered to the Earth radius. Is simply the theory 
of electromagnetic wave propagation In a horizontally stratified conducting 
aedlua, whose conductivity Is a scalar* Obviously, this statement applies 
only to Che signals propagating directly from the source Into the strato- 
sphere; the subjects of Ionospheric reflection, Schumann resonances and tensor 
conductivity In the upper atmosphere are outside of the scope of this chapter, 
the primary purpose of whicn is to discuss the stzacospherlc sferlc signal 
structure of the longest time scales (> 1 sec). Even so, the general electro- 
magnetic character of Cha theory of sferlc propagation in the stratosphere Is 
not widely appreciated. This is due primarily to the traditional development 
of Che theory, emerging from the Green's function theory of the electrostatic 
signature of charge-separated cloud structures (Holzer and Saxon, 195''.) to 
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coiKid«racioiia of tima-dapondcnc aloccrlc signaCure* of Che cloud atrucCures 

\ 

(Aoderaon and Freiar, 1969; lUingvorch, 1972). These scudles. some recognl*- 
slog Che significance of Che Green's funcclon approach, Jusclflablj slapllfled 
Che Cheory by explolclog Che coodlclon of curreac conservaclon and by Ignoring 
Faraday's lav of Induction Co ellialnaCe Che oagneclc field. Such theories are 
not sferlc theories In the proper sense because Che signal, not being an 
electromagnetic vave, is not radiated. Vlch the advent of measuring lonospha** 
rlc potential differences at balloon altitudes (Mozec and Serlln, 1969, Mozer, 
1971a and 1971b), Interests In Che cransmlccal of Ionospheric or cropospherlc 
pocenclal differences chrough Che upper atmosphere were enhanced In the early 
parts of Che past decade (Raid, 1965; Kellogg and Weed, 1969; Atkinson et al. , 
1971; Volland, 1972; Park and Oejnakarlncra, 1973; Bostrom and Fahleson, 1974; 
Chiu, 1974). These studies, some Involving Che detailed conductivity scruc~ 
Cure of Che upper atmosphere, are primarily concerned with the crans^sslon of 
sinusoidal special and/or temporal fields of given scale; therefore, only 
limited Incerprecaclon on Che source and stracospherlc signal structure of 
sf erics can be obtained from these studies. Nevertheless, the enhanced inter— 
eat sclmulaced mors considerations of sinusoidal structures (Dejnakarlntra and 
Park, 1974; Bostrom and Fableson, 1974) in which Che importance of Che elec— 
cromagneclc character of sferlcs ac altitudes below - 70 km was recognized. 
Therefore, It remains for us to bring this development Co Its next logical 
step, which is Co provide the formal and seml-o;aanCitaCive mechanism for Che 
link between source structures and the structures of stratospheric sferlc 
signals. 

Since our problem, as defined above, concerns the radiacion of the sferic 
signal in a conducting medium, sources arc an integral elemenc in the formula- 
tion in terms of Maxwell's equations. We assume the stratospheric me »iua to 
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b« aonpolarlzablt cod. ooaa«gn«eizabl«v to tuT* cha panlclvlty And panMAblll- 
cy of fr«« tpaco, tod to hmv* a horlzoncally atratlflad scalar conducelylty 
o(z). Elactrlc and aagnaclc fialds (?» ?) arlaa as cba rasult of a distrlbu- 
tlon of sourca chazgas (r, t) and sourca currants (r, t). Thasa fialds 
and sourcaa ara govamad by Maxwall's aquations: 

Causs* lavs: c7*f*Q (1) 

O 8 


7*8-0 


( 2 ) 


?araday'f law: 7 x ^ ~ - 0 (3) 

Aapara's law: 7x1 - of ^a^*o 


In (4), of la tha coodtctlon currant in tha stratosphara aasoclatad with tha 
radlatad flald f» Froa thasa Hazwall's aquations, ona darlTss tha couplai 
radiation flald aquatlona In tarns of a vactor potantlal A and a scalar poten- 
tial # such that 
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(5) 
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( 6 ) 


Tha potentials ara go^emad by tha coupled inhotaoganaoua aquations 
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Tb«M •quaclons h«T« Cbe dlMd^ancag* of bo log ^ncraccobl* «von if a convani*' 
one gauga for X Is ehosaa. Indaad, (7) and (8), whan aolwad, would Tlald ouch 
mora inforaaclon Cha currancly awailabla straioapharlc sfarlc data, slnea 

only sfarlc alaccric flald aaasuraaancs ara awallabla. Thus, wa use an azacc 
but slightly iiff treat reduction of Maxwell's aquations (P * (4) to obtain a 
singla vactor aquation for tha electric field* By talcing the tiae derivative 
of (4) and subsequent use of (1) and (3) to eliminate the magnetic field, 
arrive at the stratospheric sferic equation 
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Tha* left band side of (9) contains Cha aama linear operators as the so called 
'telegraph equation', so the f field of (9) is unquestionably Che causal 
radiation field generated by the sources cn the eight hand side of (9). The 
magnetic field, satisfying (2) and (3), is obtained from E by tiae integracion 
of (3). Our discussion of sferic propagation will be based on (9)* 

Irrespective of the specific solutions of (9), some general properties of 
the solutions can be gleaned froa an examination of tha terms in Che equa** 
Cion. The first two terms on the left bend side of (9) constitute the caut^l 
electromagnetic wave operator, while Che third term repres<*nts the effect of 
f ialte-coDductivity danplng with a damping tine 
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This dsaplag claw T, s fuaccloa of slclcud*, rspreasocs ch« aaxlasl ciaa seals 


I 

of asdiua nspoosss Co lapulslva sourcss. Uhlls s spscific dcaoascracioa of 
Chis sffccc will bs concrasced vlch obssrraclons in s lacsr tscclon, Che 
^ cypical dspandance of t and o upon Cha alcicude z is shown on Figure 9. From 

K Chls figure » we -«bserve chac chls anirt ms 1 time scale decreases wlch alcl** 

cude* Thus, if we imagine Cba propagacion upwards of a complex signal scruc- 
^ Cure, coocaining a specCJust of frequencies, Che lover Crequencies of cne 

I pacIceC are allowed ac low alcicudes while ool 7 Che higher frequenc 7 componencs 
can pcnecraca Co higher aleiCudes. In a sinusoidal wave reprasencacion, chls 
efface is of can ref acred Co as adaiccivlC 7 (Dejoakarincra and Park, 1974). 

Next, the source Cerms on Che right hand side of (9) need Co be discus- 
sed. Ic is generall 7 accepted chrc Lighcnlng discharge currencs generace 

sf erics, buc iC muse be noted also that the source Cara is not 1 but Che 

s 

clae race of ciiaoge of 5 ; therefore, a correlacive scud 7 of Che sferlc elec- 
crlc field signal structure with Che discharge current structure mist cake 
this feature into account. Another interesting feature of (9) Is that cha 
gradient of charge distributions in Che chundarr’To'id is also a source of 
sferlc radiacion. Indeed, sonr slerlc signals with poLacicy reversals appear 
to be scrlklnvlj •laliar in structure Co polarity reversals of pocentlal 
iiClencs in the chimdarscoca cloud (Vonnagut ec al., 1966). In any case, 


wlCh two possible sferlc sources Co compil.aL'a saccars, the presently avail- 
able data do noc permit an exhauaclve study of che source sCi*ucCure; chcra- 
fore, we shall acceape Co discuaa Che propagacion sciniccure by considering che 
properciea of Green's funcclon solutions Co che sferlc equation (9). 

Equation (9) is a linear scalar differencial equacion on a vector field 
with a vector source, which must be defined in cerus of three independent 
scalar components functions although che three components are not coupled. 


C.Zu 






Th.la iaportanC f«acucc is espaci&lly imporcaat in ths consideracioa of bounds* 
ry coaditioas for E. For axanple, in Che limit of a parfectly coaduccin«; 


earchf the Green* i fimction for the component in a cylindrical coordinate 


syatcm (p, z, oust satisfy the boundary condition ■ 0 at the Earth's 


surface z - 0. Since the OMthods of constructiixg Green's functions satisfying 
such boundary conditions are standard [see for example Bolzer and Saxon (1952) 
for the image source method, and Chiu and Hilton (1977) for integral equation 
methods], and since ve are at present uncertain of the stnicture of the vector 
sources, we shall deal with the simplest elements of sferic propagation 
Green’s function — that of infinite boundaries. 


In view of the above discussion, it suffices for us to investigate the 


properties of the Green's function G ( r, r^, t, f^) fcr the scalar equation 
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Th« cooductivtty 0 ( 2 ), with a typical altitude dependence shown on Figure 9, 
can be parameterized (below ^ 70 ka) aa 


0 ( 2 ) ■ a exp (z/h) 
o 


( 12 ) 


where o “ 5 x 10”^^ mho/m and h » 6 taa. The Green's function G for (11) and 
o 


(12) has not been considered although, as mentioned above, the siausoidal wave 
solutions of (11) have been considered by many authors. Ue consider this 
problem In the limit of azimuthal symmetry (no ♦ dependence) in a cylindrical 
coordinate system (o, z, >). 
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Using chs delta function representations In sinusoidal and Bessel function 
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• -l(i»(t-t ) 


o du 


7 «(x-x^) - // Tdy JpCsy) JpC^^y). 


and cransforoing Cha coordlnaCa z Co a nav coordlnaCa Z = axp (z/2h) we can 
%rrlca an exact caprcsancaclon o£ G by Cha dacompoalCional aachod. [saa for 
axaapla Chiu and Hilcoc. (1977)] : 


G r *10110 I *\dA / ^ J Ocp)J (kp )J (X^)J (X^ )/ 

r ^o -'o ^ o o o p P p 


(\^HL4u a h""ui) 
0 0 


2 2 2 2 2 

where p"" • 4h (k • uiVc )• Because (16) involves integration over the order 
of Bessel functions , it is not tractable for the purposes of data Intsrpreta* 
tlon even though the X-integration of (16) can be carried out exactly by 
contour Mthoda [see Watson (1966) p* 429] a 

^^e salient properties of (11) and (13) nay however be discuased by other 
nathods involving Green's function solutions* One such taethod is to consider 
a simplified but tractable version of (11) and (13) while keeping la mind that 
the simplification must not affect the essential properties of the solution. 
Such methods of considering a *surrogate'* have been applied in many areas [see 
for example Chiu (1976)]. In the case of (13), we note that the character of 
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cb* •quaclon involve* cauMl propegetloa of a wave signal as well as “diffu- 
sion’ damping from Che Linear time derivative term, which is Che source of 
solution diffictilcies. To obtain a “surrogate" we note that if the conducti- 
vity is constant the causal propagation of the sferlc wave is unaffected but 
we may insight to the sferic signal stnicture by eliminating mathematical 

difficulties. 

If we assume the conductivity to be constant, an exact solution to (11) 
and (13) ran be written in closed form (Morse and Feshbach, 1953, p. 868) 


,(- Vj) . c ) - RI t 

o 


+ J. [(Mj) Wq /r^ - c^(t-t^)"] 9 [c(t-t^)-R]} (17) 

2 V - c^(t-c )^ 
o 

2^2 

where R = (p-p )* -i- (z-t ) , and 9(x) • 1 if x > 1 and zero otherwise- 

0 . O ' 

The character of Che Green's function (17) becomes apparent if we consi- 
der Che physical basis of each of the terms. The term 5 [c(t-tQ)-R]/R is 
simply the Lienard-Uiecherc potential of the causal electromagnetic radiation 
field. This is expected if sf erics, at whatever frequency, are radiated 

signals. The finite conductivity of Che medium gives rise to two ef- 

fects. First, the entire signal structure is damped in time by Che factor exp 

1 2 2 
[“ ( V 2 )u o c" (c-c.)]. Thus, T i 1/u 0 c“ sets the upper limit to the time 

0 o o 

scale inspite of the source time scale implied in the Lieoard-Weicherc term. 
This effect 1s expected whether T is constant or rot, and it has been dis- 
ciissed previously la this section. In addition, Che "diffusive" effect of the 
conduction current Introduces a second, tall-liite element in the signal struc- 
ture in Che form of the Bessel function term of (17). It should be noted that 
this effect is a tall-like signal because Che causal condition (signified by 
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ch« d-functloo) raqulres ioagloary arguneacs for th« 3«s el fuacCioa* Th* 
tlo« seal* of chls tall Is qoc exactly t but Is of the same order (~ a few 
seconds in the stratosphere), as we shall numerically illustrate in a later 
section. This tall-lllce structure may be a very important feature in data 
interpretation because, In contrast to madia wave responses to Impulsive 
disturbances, the sferlc response. In the stratosphere consists of a radiation 
pulse plus a single tall. This structure is boms out by observations In the 
preceedlng section and discussed In more detail in the next section. 

Next, let us Investigate another "surrogate" of (13) which emphasizes the 
vertical structure, as Imposed by o ■ exp (t/h) , but deemphaslzes the 
causal element of the signal structure. In (17), the ciusal functions, d 
[c (t-tg) - a] and 0 (c(t-tg) - al, are formally very Important, but for 
practical Interpretation of signals received ■ 100 km from the source the 
causal time delay is only of order of milliseconds or less. Therefore, for an 
instruffleot with response time longer than milliseconds it may be profitable to 
consider the Impllcstlons of the vertical conductivity structure at the ex- 
peusa of maiataining causality. This means that we consider the Green's 
function In the "diffusion" limit of (13), l.e., c ■*■ and the second time 
derivative, term of (13) vanishes. In this limit, we obtain a Green's function 
representation 
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carried out 
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with the Identity 


(z/2b). 


The p— integration in (18) can be 
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Cha Gr««n's funccioa (18) yield* Che claeelcel 'diffusive* elemeacs of Che 
fore 
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Noce c^ec ■ o(z). In Clie screcoephere where » 9 X 10**^2 aec"l 

# 

C 10, wc.aocc chac Che arguaenc of l2hit ^ (20) Is very small for ■call- 
like" clae ecales c-c^ - 10 sec. Thus, for Che "call" pare of Che signal 
scrucCure, we can use Che approxlmaclon 


^2hk " (x/2)^ ■ exp [2hk In (x/2)] . (21) 


Applying (21) Co (20) and noclng chac Che argximanc of I-^^^ Is less Chan unicy, 
we obtain an approximace expression for Che "tall" part of G by carrying out 
Che Integration In k in (20) 
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where 


a “ 2h|in[a^CC /4(c-t )] I “ r 
0 0 


( 23 ) 
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Comparisoa of (22) Co tha "tall" part of (17) shows a number of revealing 

features despite the loss of tne propagation parameter c in (22). First, the 

causal factor S(t-Cg)/(c-tQ) in (22) directly corresponds to the causal factor 

d(c(c-t ) - R]//(P.^“C^(t-g)l in (17) in the limit c ♦ The damping factor 

2 2 

in (22) has a time scale which is proportional to l/o(z) (because a Z 

CC (z)) In (17), although because of the loss of c in (22) the rwo time 

scales are very different. Nevertheless, we regard this feature as indication 

2 

chat Che time scale T s 1/u oc can be evaluated at various values of z even 

o 

chough (17) is derived with constant z. The scaling in the spatial structures 

of Che signals In (17) and (22) are, however, very different. With Che use of 

(23), we note that the spatial struccure is roughly proportional to 
2 2 3/2 

z/[z > ()J*n) ) 1 , as is usually expected in Green's functions of the 

o 

operator 7^. The spatial structure of (17), for constant a, is roughly 
1/R. However, if one were Co evaluate (17) with a~o(z) Chen Che magnitude of 
Che "tail" structure will be dominated at high altitudes by the o(z) factor in 
the numerator. Considerations of (22) shows conclusively chat this i? <jnphys- 
ical. In ocher words, the evaluation of (17) with the artificial substitution 
o*o(z) is roughly justified in determining time scales and horizontal spatial 
scales of the signal structure but is not justified in the incerpretation of 
Che magnitude of response as fxinction of altitude. The altitude dependence 
t/lz" + (ij-o^)"! may be more truly representative of Che signal structure. 
In this section, we have accenpced to review what has been done in Che 
area of sferic signal propagation in Che stratosphere. have, in Che pro~ 

cess, carried the analysis one loglcsi seep further by considering the problem 
in terms of radiation and propagstion in a conducting medium. The intimate 
lacerpiav between the space-time structures of the source and of the propaga- 
tion medijo is shown to be inCerpretable Ln terms of lireea's functions, ’-hile 




it is rather unfortunate that an exact representation of the aferlc Green's 
function requires numerical work far beyond what is appropriate for this book, 
investigation of some 'surrogates'’ reveals properties which are borne out by 
balloon observations of sferic electric fields in the stratosphere. 
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DISCUSSION OF TIME CONSTANTS AND CURRENTS 


In ch« Last few aecCions w« have preeeaced a zoo of sferic observacioos 
in Che scracosphere and Che siachemeclcal basis co help underscand cheir propa~ 
gatlon cheracceristics. Now ue will use Chls cbeory anu Chese observacioos in 
a few quancicacive calculaCions in which we will accempc Co explain some of 
Che sferic signacure characcariscics and Co calcolace currenc densicies asso- 
ciaced wich sferics. 

The cxacC Creen's fimccion given in (17) has been approximaced in (22) by 
assuming conscanc o* Equacion (22) was Chen evaluaced ar. z*30 ^ using o as 
given by equaCion (12). The upper plot in Figure 10 shows daca from 0016 UT 
August 14, 1974 which included a small sferic chac did noC sacuraca our elec- 
cronies (also see Figure 4). On cop of Che data is Che ploC of (22) after 
normalizacicn of Che ampliCude Co 3.54 voles at c^. Since (22) assumes a 
conscanc a (evaluated ac z*30 km) ic cannot be used Co directly calculate Che 
ampliCude of Che expected signal in Che scraCosphera as discussed in the I ist 
section. However (22) does give a good qualiCati'/e explanation of the sferic 
cail. The sferic Calls seen in *1 1 of Chese data are shorter than the re- 
sponse else of Che medium as shotm in Che bottom panel of Figure 10. This is 
a conductivity aeasuremenc on Che same flight as Che upper data taken about 3 
hours later. In this aeasuremenc che upper vertical probe is incemally 
shorted to payload ground for four seconds (c ■ “4 to 0 in Figure 10) and Chen 
refloated. The time response of che medium is inversely related co che ambi- 
ent conductivlcy a and in this case an exponential fit to che data after 
refloating gl/es a conductivity a • 6.2 x lO”^" aho, meter with probe effects 
incluaad. This incidentally compares -very favorably to (12) which, when 
evaluated ac z-30 'cm, gives j(30 km) ■ 7 .j* i iiho/ a. For a more dccal^ed 


vohs/meier 
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seconds 

KlV,iiro MK (rppt*r panel) About 10 seconds of raw rlectrlc ii.it. < — 

surrounding, the sferlc observed at OOih.S IT, Aupust 1-*, 1^74 (m*.* 

also FI pure ^), Plotted over the data In the upper p.inel I the 
r.reen's function fron equation (22). The hot ton panel Is \ con- 
duct I vl tv measu»*enent sht^wlnp time responsi* of the mediun t«» be 
lonper than the tall on the sferli hv about a second. See text tor 



discus* Ion of this type of cooductlvicy aetsurement tn the stratosphere see 
Holzvorth (1980). 


* ^ ' 


To emphasize the fact that the sferic tail really does decay faster than 


the ambient time response, the upper and lower data were both fit with expo- 


oentlals resulting in time conatants of i * 1*56 seconds for the sferic 


(upper) and 2.86 sec for the conductivity measuxements (lower). This shows 


that the sferic signal Is indeed not simply an ohmic current but that the 


mathemacical development in Section IV really is necessary to explain the 


propagation effects. It should also be noted that the final sferic return 


tails (from negative value near -4 volta track up to base line values of 


Figure 4 through 7 seen in the horizontal probes (lower panels with and 


are due to local time reiponse effects after the ground fluctuation. The time 


constants determined from and follo%ilD^ large sferics agree quite well 


w1.th our probe shorting time responses. 


The actual sferic current density flowing past the balloon cannot be 


unraveled easily from the electrical signature because of its nonob'-iic propa- 


gation effects. However^ dramatic and complex ground fluctuations seen in /til 


large sferics in Figures 4*7 and also by Bering (private conmninication, 


1979) do result in redistribution of stratospheric ions and hence currents. 


In most cases the horizontal currents in the sferics are tne same order of 


magnitude as the vertical. 


The observed sferic signatures can be used, how-^ver , to place an upper 


Unit on the current density and total current flowing in the stratosphere 


during a sferic. If the entire transient of the vertical field is assumed to 


drive an ohmic current it would have an average peak current density of J. 


■ 6.2 X 10*^^ mho/m * 5 v m - 31 pico anps/n. Thus 31 pico aaps n" 
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Is an uppsr liaic to tha pask sfarlc currant dausicy which alghc ba oucli 
smallar. Coaparlng this uppar Limit to tha actual dc currant aeasured over a 
chuiid€rstora* which cm b« av€r 40 plco amps/a^ for hours (se« Figure 2 and 
HoLzworth (1980)), tha data show that afarlcs arc probably not vary important 
In charging up tha lonoaphara ralatlva to tha ground but that tha dc dynamo 
cxirrants driven by thunderstorm convection moves anich more total charge than 
all tha sferlcs together. 

This chapter has provided an observational and theoretical background for 
tha study of stratospheric sferic signals in tha frequency range below 2 He. 
We have described the balloon borne high impedance electronic instruments used 
to detect tha electric signals in this range and presented summary sferic data 
from dozens of balloon flights over and near thunderstorms. We then extended 
tha theoretical background to treat radiation of sferic signals generated in 
the lower atmosphere with propagation into the stratosphere. This mathe-* 
matical derivation of an electrical sferic equation resulted in an exact 
equation which was then simplified for sample calculations in order to help 
expLaii the temporal structure of the observed sferic signals. The observe* 
tions in Section III can not be used as sferic source Locations without a 
leccnd balloon measurement so that the amplitude attenuation with distance can 
be calculated. 

In this '^sL section we applied the Green's function theory/ developed in 
Section 17 in a sample calculation which helps to explain the time constants 
in the observed sferics. An upper estimate for the peak current density in a 
sferic was calculated assuming the observed sferic electric signature was all 
due to ohmic current. This is likely to be a gross overcstlaare but neverthe- 
less shows that Che total current carried to the ionosphere as a result of 
sferlcs is probably iaslgnif leant in comparison with the upward dc current 
over thunderstorms . 


( . ^ 


Much <*pn be gainad in a further study of sfarlca by higher ciaie rasolu** 
tioa inatruiaentacloa and eulciple slaultanaoua observationa. Also note accu- 
rate approxixaaclona to equation (17) vould allow more useful comparison to Che 
observed sf erics, fvfcrlc resea.rch in Che scratosphere has JusC begun and this 
chapeer atceopcs Co tie soma of Chis early work cogecher. 
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Abstract 


Based on the occurrence frequency of wave-like structures in some 24,000 
measurements of thermospheric density, we demonstrate that there is no signi- 
ficant hemispherical correlation between satellite-measured thermospheric 
wave-like structures and the tropospheric vorticity area index in the epoch 
1974-1976. However, in the northern polar zone, (60“-90*) N, wave occurrence 
frequency is greater for large values of the 500 mb vorticity area index. 
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introduction 


In two previous studies of the geographical distributions of wave^llke 
density disturbances In the thermosphere, It was established that the dlstrl‘~ 
butions show patterns of occurrence which match persistent features of the 
troposphere, such as regions of strong wind shear associated with the subtro-* 
plcal jet stream (Rice and Sharp, 1977) and seasonably persistent regions of 
deep convective activity Indicated by high occurrence rates of lightning at 
dawn (Chiu et al«, 1979). Such studies, together with studies of the auroral 
origin of large-scale thermospheric wave-like density disturbances (Potter et 
al«, 1976), offer suggestions as to the . origin of persistent global-scale 
sources of thermospheric gravity waves, which we identify a priori with the 
thermospheric wave-like density disturbances discussed in the above studies* 
At the small-scale end of the spectrum of atmospheric disturbances, thermo- 
spheric waves associated with local-scale (mesoscale) disturbances of the 
troposphere, such as thunderstorms, tornadoes and hurricanes, have also been 
reported (Kung et al., 1978). Thus, considering the entire spectrum of atmo- 
spheric disturbance scales as a whole, one is naturally prompted to ask if a 
similar correlation between satellite-measured thermospheric disturbances and 
transient synoptic scale (weather) systems in the troposphere cay not also 
exist. 

The existence or non-existence of such a correlation has significant 
implications for the viability of some suggested physical mechanisms for the 
so-called "sun-waather” effect (cf. Wilcox fet al., 1976); in particular, those 
mechanisms requiring transmission or modification of an atmospheric disturb- 
ance between the upper and lower atmosphere will be affected. For example, 
Hines and Halevy (1975) suggested that a “solar influence must somehow modu- 


late the meteorological noise" to produce the "sun-nfeather" signal. A condl~ 
tlon for such physical mechanisms to be viable is that there exists a rela** 
tlonshlp between "meteorological noise" of the proper scale in the upper and 
lower atmospheres. In this correlative study of upper and lower atmospheric 
data of "meteorological noise,” we do not address the question of the "sun- 
weather" effect; rather, we focus upon the correlation, if any, of satellite- 
measured upper atmospheric "meteorological noise" with that of the troposphere 
as is indicated by the vortlclty area index (VAX). Although we demonstrate 
that there is no significant hemispherical correlation between satellite- 
measured thermospheric wave-like structures and the tropospheric vorticity 
area Index in the epoch 1974-1976, we cannot make any comments on the "sun- 
weather" effect in general because it is possible that the solar signal can be 
carried by upper atmospheric "noise" in portions of the scale spectrum not 
covered by our satellite data. In the lilgh- latitude region (60®-90“) N, 
thermospheric wave occurrence does depend upon the vortlclty area index. If 
this correlation is to be identified with the "sun-weather" effect, then the 
likely mechanism would be thermospheric processes associated with polar pheno- 
menon such as the aurora. 
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DATA AND ANALYSIS 


The correlation study Involves two sets of data. The lower atmospheric 
data set Is the vortlclty area Index (VAI) for 1974-1976 compiled by R. G« 
Olson, W, 0. Roberts and E. Gerety (Shapley and Kroehl, 1977). The VAI 
values. In units of 10^ km^, are derived from the 500 mb height analysis 
records of the U. S. National Meteorology Center for the northern hemisphere 
above 10° N, and for both 0000 UT and 1200 UT. The upper atmospheric data set 
consists of 24,000 density measurements In the northern hemisphere some of 
which contain wave-llke structures. These data were obtained from cold 
cathode Ion gauges onboard Atmospheric Explorer spacecrafts A£-C, AE-D, AE-E 
and the Air Force spacecraft S3-1. Data obtained In the altitude range 135- 
225 km are used, k description of the cold cathode Ion gauge Is given In Rice 
et al. (1973). Figure 1 of Reference A (Chiu et al., 1979) depicts a rela- 
tively disturbed neutral density profile as a function of altitude; similar 
wave-llke structures observed by other Instrvuaents onboard AE-C have been 
interpreted as gravity waves (Reber et al., 1975). In this figure, four 
classes of wave-llke structures are illustrated; the four wave classifications 
are identified with an index value running from 1 to 4 (Rice and Sharp, 
1977), The classifications are: I ■ 1, indicating a smoothly varying pro- 
file; I ■ 2, indicating minor activity with up to 5% density amplitude, I ■ 3, 
indicating waves of 5-15% density amplitude peak to peak; and I ■ 4, indi- 
cating waves of density amplitude greater than 15%. For this study, wave 
classifications 3 and 4 have been combined as an indication of the unambiguous 
presence of a gravity wave. Although in general wave classifications 3 and 4 
have different wavelengths, we have not yet found any indication that they 
originate from different sources. For the data from AE-D, we have excluded 
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periods when the spacecraft was undergoing nutation. In addition to data used 
in Reference A, data from the AE-E satellite have been added to the data base, 
comprising a total of 23,691 density measurements at 5 km altitude Intervals 
of the above satellites in the northern hemisphere in the period 1974-1976. 

The thermospheric index is the wave occurrence frequency (WOF) which is 
defined as the percentage of wave measurements of classes 3 and 4 found in the 
portion of the total number of measurements (23,691) which corresponds to the 
stated conditions of the particular correlation parameter. Thus, if a certain 
WOF is 0.24 for a stated correlation condition, then classes 3 and 4 waves are 
found in 24% of the density measurements corresponding to that condition. 
Generally, if a relationship exists, the WOF will depend on the VAI; and this 
relationship ’s quantified in a plot of WOF as function of VAI. The error 
analysis for this relationship is particularly simple if we assume normal 
distributions. If waves are found in N samples of a total of M appropriate 
measurements, the WOF is (N/M) ± (/N/M). The results of the correlations are 
stated in WOF because satellite measurements are not evenly distributed in 
space and time and VAI mean values also vary secularly in the period 1974- 
1976. 

1. Hemispherical relationship between WOF and VAI. Figure 1 shows the 
relationship between WOF and VAI for the northern hemisphere for three time 
shifts between the indices. The numbers near the bottom of the figure indi- 
cate the distribution of measurements in each VAI Interval for 0 day time 
shift (measurement distributions for other time shifts are similar). Evident” 
ly, VAI values in the sampling period do not form a normal distribution. 
Nevertheless, the significant result of this figure is chat thermospheric wave 
occurrence frequency (WOF) in the northern hemisphere is * (25 ± 5)%, indepen- 
dent of VAI values. Even though the 5% fluctuations shown are outside of 
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WAVE OCCUR 



Kij^ure 1. Correlation between wave occurrence frequency in the northern 
hemisphere and the vorticity area index* 
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nortaally distributed statistical limits, possible dependence of the WOF upon 
the VAI at the < 5% level is hardly of significance. 

2. Zonal relationship between WOF and hemispherical VAl. Figure 2 

shows the relationship between WOF in the polar (60*N-90“N, P), mid-latitude 

(30*N-60“N, ll) and equatorial (0*N-30*N, E) zones and the hemispherical VAI at 

5 2 

0 day time shift. Except for statistical fluctuation in the (40-50) x 10 km 
bin of VAI, the WOF in the equatorial (E) and mid- latitude (M) zones (and/or 
the sum of the two) is again independent of VAI, at a value of “ (22 * 5)Z. 
The WOF for the polar (P) zone shows a significant (factor of two) jump 
from - 15S: to - 303: for VAI values greater than 30 x 10^ Similar depend- 

ence is iound for ± 2 day tine shifts. Examination of the measurement distri- 
butions among the three zones (bottom of Figure 2) indicates that the centroid 
of VAI for the polar zone (- 50 x 10^ km^) is higher than that of the other 
two zones (- 25 x 10^ km^). This prompts us to note that most polar zone data 
come from the polar— orbiting AE— D in the epoch (Nov. 1974 - Jan. 1975) when 
VAI values happened to be high because of the northern winter. Thus, the 
strong dependence of the polar zone WOF upon the VAI may be due to preferen- 
tial sampling of wintertime VAI by the AE-D satellite at high latitudes; 
nevertheless, it cannot be denied that WOF is also unusually high during this 
interval Irrespective of sampling. 

Because of the abnormal measurement distribution, we examine this “lati- 
tude” effect further by contrasting the AE— D data with the AE-C data, which 
also contribute to the polar zone. The AE-C data extends up to only 68® N. 
Figure 3 shows the polar zone WOF for the AE-C and AE-D data as functions of 
VAI at zero time shift. The centroid of the measurement distribution for the 
AE-C data set is at VAI » 35 x 10^ km“ while that for the AE-D data set is at 
VAI » 50 X 10^ km^. Both data sets show a substantial increase in the polar 
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WAVE OCCURRENCE FREQUENCY IWOH IN LATITUDE ZONES 
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VAI, 10® Km^ 


Figure 2* Correlation between wave occurrence frequency in the equatorial, 
mid-latitude and polar zones with tha vorticity area index* 
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WAVE OCCURRENCE FREQUENCY 



3* Correlation between wave occurrence frequency In tlu* polar 7.otu? 
with tb,e vorticlty area index for the AE-C and AE-D data sets. 
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zone WOF as VAX Increases past 30 x 10^ Therefore, it Is likely Chat the 
polar zonal relationship between WOF and heaispherical VAl shown on Figure 2 
is really a latitude effect and not due to sampling bias of the AE-D satel~ 
llte* Even if we discount Figure 3, we emphasize that the polar zone data 
clearly shows dependence of WOF upon the VAX, although one would not then be 
able to say whether it is due to a latitude effect alone or due to a secular 
increase of the average VAX during the AE~D sampling period* 

Xf this polar zone effect has some relationship with the "sun'-weather" 
effect, then it would argue for upper^lower atmospheric wave coupling assocl'- 
ated with polar phenomena such as auroral processes* We offer this as a 
speculative comment since the "sun^weather” effect need not be germane to our 
specific discussion of upper and lower atmospheric wave correlations* 

3* Hemispherical relationships between WOF, VAX and interplanetary 
magnetic sector boundary crossings* In Figure A, we show the northern hemi~ 
sphere WOF and VAX, averaged over the measurement distribution (shown at the 
bottom of the figure), within ± 6 days of interplanetary magnetic sector 
boundary crossings in epoch 1974-1976* Even though the average VAX is taken 
over an uneven measurement distribution, tiie characteristic signature of VAX 
dependence on days from sector boundary (Wilcox et al., 1976) is found. The 
hemispherical WOF, however, remains at *” (25 ± 5)X, Indepe idc Jt of days from 
sector boundary crossings* We have also examined the relationship between 
polar <;one WOF and the interplanetary magnetic sector boundary crossings* 
Again, we found that the polar zone WOF to be independent of days from sector 
boundary crossings, although we have only a small data aet to work with* 


D-13 


AVERAGE 




CONCLUSIONS 


Based on the occurrence frequency of thermospheric wave structures In 
some 24,000 measurements of density, we demonstrate that there Is no slgnlfl* 
cant hemispherical correlation between thermospheric waves and the tropo* 
spheric vortlclty area Index In the epoch 1974-1976. Further, wave occurrence 
frequency does not exhibit the characteristic feature of the VAI In relation 
to Interplanetary magnetic sector boundary crossings. In the polar zone, (60* 
- 90”) N, wave occurrence frequency does Increase significantly if the heml- 
spherical VAI Increases above 30 x 10*^ km . In contrast to previous work 
(Rice and Sharp, 1977; Chiu et al., 1979), which correlated thermospheric wave 
occurrence to geographically and seasonally persistent features of the tropo** 
sphere, this work focuses upon the study of transient (weather) systems as 
represented by the VAI. Consequenly, the null relationship between VAI and 
the wave occurrence frequency outside the polar zone may reflect the dominance 
of the persistent sources In the non-polar troposphere. 
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